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ABSTRACT 

We present an in-depth study of the hard-binary population of the old (7 Gyr) open cluster NGC 
188. Utilizing 85 spectroscopic binary orbits out of a complete sample of 129 detected binary members, 
we study the cluster binary frequency and the distributions of binary orbital elements amongst the 
main-sequence, giant and blue straggler populations. The results are derived from our ongoing radial- 
velocity survey of the cluster, which spans in magnitude from the brightest stars in the cluster to 
V = 16.5 (about 1.1 - 0.9 M Q ), and extends to a projected radius of 17 pc (~13 core radii). Our 
detectable binaries have periods ranging from a few days to of order 10 4 days, and thus are hard 
binaries that dynamically power the cluster. The main-sequence solar-type hard binaries in NGC 
188 are nearly indistinguishable from similar binaries in the Galactic field. We observe a global 
solar-type main-sequence hard-binary frequency in NGC 188 of 23 ± 2 %, which when corrected for 
incompleteness results in a frequency of 29 ± 3 % for binaries with periods less than 10 4 days. For 
main-sequence hard binaries in the cluster we observe a log-period distribution that rises towards our 
detection limit, a roughly Gaussian eccentricity distribution centered on e — 0.35 (for binaries with 
periods longer than the circularization period), and a secondary- mass distribution that rises towards 
lower-mass companions. Importantly, the NGC 188 blue straggler binaries show significantly different 
characteristics than the solar-type main-sequence binaries in NGC 188. We observe a blue straggler 
hard-binary frequency of 76 ± 19 %, three times that of the main sequence. The excess of this binary 
frequency over the normal main-sequence binary frequency is valid at the >99% confidence level. 
Furthermore, the blue straggler binary eccentricity - log period distribution is distinct from that of 
the main-sequence at the 99% confidence level, with the majority of the blue straggler binaries having 
periods of order 1000 days and lower eccentricities. The secondary-mass distribution for these long- 
period blue straggler binaries is narrow and peaked with a mean value of about 0.5 M Q . Predictions for 
mass-transfer products are most closely consistent with the binary properties of these NGC 188 blue 
stragglers, which comprise two-thirds of the blue straggler population. Additio nally we compare the 
NGC 188 binaries to those evolved within the sophisticated iHurlev et ail (|2005| ) iV-body open cluster 
simulation. The main-sequence hard-binary population predicted by the simulation is significantly 
different from the main-sequence hard-binary population observed in NGC 188, in frequency and 
distributions of period and eccentricity. Many of these differences result from the adopted initial 
binary population, while others reflect on the physics used in the simulation (e.g., tidal circularization). 
Additional simulations with initial conditions that are better motivated by observations are necessary 
to properly investigate the dynamical evolution of a rich binary population in open clusters like NGC 
188. 

(galaxy:) open clusters and associations: individual (NGC 188) - (stars:) binaries: spectroscopic - 
(stars:) blue stragglers 



1. INTRODUCTION 

Binaries influence the evolution of open clusters and 
are catalysts for the formation of anomalous stars like 
blue stragglers (BSs). Stellar dynamicists have long rec- 
ognized that encounters involving binary stars can supply 
the inevitable energy flow out of the cores of star clus- 
ters. These complex dynamical dances also lead to stellar 
exchanges, binary orbit evolution, close stellar passages, 
mass transfer, mergers, and stellar collisions, all of which 
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lead to an array of new stellar evolution paths and prod- 
ucts. Observations of the old (7 Gy r; ISaraiedini et al.l 
119991 and see also lMeibom et al.ll2009[) open cluster NGC 
188 reveal a rich binary population and a wide variety 
of such interaction-product candida tes, including BSs , 
sub-subgiants and X-ray sources (e^g JBelloni et al.ill998l : 
lGondoidl2lMlGeller et aIE(M. l2009h . This abundance 
of anomalous stars may be a direct result of the dynam- 
ical evolution of the single and binary stars as they in- 
teracted over the lifetime of the cluster. 

The field binaries are an important comparison popula- 
tion to open cluster binaries, as we would not expect field 
binaries to undergo dynamical encounters due to the low 
stellar density of the Galactic field. Therefore compar- 
isons between open cluster and field binary populations 
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may provide insights into the impact of stellar dynam- 
ics on shaping a binary population and creating exotic 
stars. There have been a number of studies of the field 
dwarf binaries in the past (notably, lAbt fc Levvlll976l 
iDuquennov fc Mavorlll99ll ) . Currently the most compre- 
he nsive and comple t e sury ey of field dwarf binaries is that 
of iRaghavan et all (|2010l hereafter RIO). They analyze 
the multiplicity of a complete sample of 454 ^F6 - K3 
dwarf and subdwarf stars within 25 pc of the Sun, find- 
ing a frequency of multiple systems of 46% (34% ± 2% 
strictly binaries). The orbital period distribution is fit 
by a log-normal curve centered on log(P [days] ) = 5.03. 
For binaries with periods longer than the circularization 
period, RIO find that the orbital eccentricity distribution 
is roughly flat out to e ~ 0.6. The mass-ratio distribu- 
tion for binaries is found to be roughly uniform between 
0.2 and 0.95, with an additional peak at a mass ratio of 
u nity (commonly ter med "twins"). 

iGeller et al.l (j2009L hereafter Paper 2) present 85 solar- 
type binary orbits in NGC 188 derived from our ongoing 
radial-velocity (RV) survey of the cluster. These bina- 
ries are of similar spectral type to those studied by RIO, 
and thus provide a comparison between field binaries and 
those that have evolved within the dynamical environ- 
ment of an open cluster. This work is an integral part o f 
the WIYN 1 Open Cluster Study fWOCS: lMathieull2000ft . 
Our observations cover a complete sample of stars within 
a magnitude range of 10.8 < V < 16.5 (1.1 - 0.9 M Q ) and 
extending to 30 arcmin (17 pc in projection or roughly 
13 core radii) 2 in radius, esse ntially out to the tid al ra- 
dius of the cluster (21 ± 4 pc; lBonatto et aLl l2005). Our 
magnitude limits include solar- type main-sequence (MS) 
stars, subgiants, giants, and BSs. For most binaries, our 
observations span ~11 years, while for some binaries our 
time baseline covers up to 35 years. Our detectable bi- 
naries have periods ranging from a few days to of order 
10 4 days. A full description of our stellar sample, ob- 
servations, data reduction routine, and precision, as well 
as a summary table listing relevant i nformation on each 
star in our sample, can be found in IGeller et al.l (2008) 
(hereafter, Paper 1). The binary orbital parameters, in- 
cluding estimates for the primary masses, can be found 
in Paper 2. 

In this third paper in the series, we study the charac- 
teristics of the NGC 188 binary population. We first 
perform a detailed analysis of our observational com- 
pleteness (Section [2]), and use the results from this sec- 
tion to study the global hard-binary frequency (Sec- 
tion [3]) . In Section [4] we discuss the eccentricity - log 
period (e — log(P)) diagram. We then analyze the hard- 
binary distributions in eccentricity (Section [5]), period 
(Section[5]), mass ratio and secondary mass (for MS bina- 
ries in Section [7] and BS binaries in Section [SJ. Through- 
out the paper, we compare the MS, giant and BS bi- 
nary populations. Additionally, we compare the NGC 
188 binary population to similar binaries in the Galactic 
field. In Section [HI we also compar e these result s to th e 
A^-body open cluster simulation of iHurlev et al.l (2005), 

1 The WIYN Observatory is a joint facility of the University of 
Wisconsin - Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 

2 We adopt a core radius of 1.3 pc (Bonatto et al. 2005) at a 
distance of 1.9 kpc, which corresponds to 2.35 arcminutes on the 
sky. 
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and we build on th e work of lMathieu fc Gel ler (2009) and 
IGeller fc Mathieul (|2011l ) to discuss the possible origins 
of the NGC 188 BS population. Finally, in Section ITuI we 
provide a brief summary and conclusions. 

2. COMPLETENESS IN BINARY DETECTION AND 
ORBITAL ELEMENTS 

There are 630 stars wit h non-zero proper-m otion mem- 
bership probabilities from lPlatais et ail ()2003f ) within the 
magnitude and spatial domain that define our complete 
stellar sample. Updating the results of Paper 1, we have 
obtained >3 RV measurements for 623 of these stars 3 . 
Thus, we have derived RV membership probabilities and 
characterized the RV variability of 99% of this complete 
sample. Two of the seven stars that don't have three 
observations are too blue to derive reliable RVs with our 
current observing setup and two are rotating too rapidly 
(both are W UMa's). The stars with nonzero proper- 
motion and RV memberships define the stellar sample 
used in the subsequent analysis, totaling 487 primary 
stars. Here, we characterize our completeness in binary 
detection and orbital solutions within this sample. 

We use a Monte Carlo approach to produce simulated 
observations of a set of simulated binaries with charac- 
teristics of the Galactic field binary population of R10. 
Specifically we choose the binary periods from a log- 
normal distribution, centered on log(F [days] ) = 5.03 
with ai gp — 2.28. We fit a wide Gaussian to the ob- 
served eccentricity distribution in R10 that is centered 
on e = 0.39 with a e = 0.31 {xl e d = 0.78) and choose the 
eccentricities from this distribution. (A flat eccentricity 
distribution fit to the R10 data results in x re d = 

2.28.) 

We use 1 M© MS primary stars and choose mass ratios 
(q) from a flat distribution such that the secondary stars 
have masses between 0.08 - 1 M Q . We do not include 
the peak in the observed field distribution at mass ra- 
tios of unity as we detect such binaries as double-lined 
systems in our spectra and identify them as binaries di- 
rectly. Once detected, we observe double-lined systems 
on a more frequent basis than the majority of the bina- 
ries in our survey (see Paper 1). The orbital inclinations 
and phases of the binaries are chosen randomly. To each 
of the simulated observations of a given simulated binary, 
we add a random error generated from a Gaussian dis- 
tribution centered on zero and with a standard deviation 
equal to our single-measurement precision of 0.4 km s . 

As defined in Paper 1, an observed star is considered 
to be in a binary if the standard deviation of its RV mea- 
surements is greater than four times our precision of 0.4 
km s -1 (or e/i > 4). In our survey, we initially followed 
iMathieul (fl983h in requiring three observations over the 
course of at least one year before attempting to catego- 
rize a star as being either single or binary. Here, we use 
our Monte Carlo method to re-evaluate this procedure. 
We determine our binary detection completeness by an- 
alyzing simulated observations of the simulated binary 
population covering time intervals defined by our actual 
observations of NGC 188. 

3 We require at least three observations before we categorize a 
given star as either single or binary. We use the term "single" to 
identify stars with no significant RV variation; a star is termed 
single if the standard deviation of its RV measurements is less 
than four times our precision. Certainly some of these stars are 
also binaries, although generally with longer periods and/or lower 
total mass than the binaries identified in this study. 
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Fig. 1. — Observational completeness curves for NGC 188 binaries resulting from our Monte Carlo analysis, 
as a function of period (left), eccentricity (P c i rc < P < 3000 d; center) and mass ratio (P< 3000 d; right), 
completeness in binary detection with the solid line and in orbital solutions with the dashed line. 
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Specifically, we have observed the stars in our NGC 
188 WIYN sample over the course of about fourteen 
years. The distribution of the time separation between 
our first two observations is bimodal with two equally 
strong peaks, one at about two months and another at 
about nine months. The distribution of the time sep- 
aration between our first and third observations is also 
bimodal, with a strong peak at about eight months and a 
smaller peak at about 2.5 years. These distributions re- 
flect the approximate time between individual observing 
runs and the number of stars that we can observe during 
each observing run (see Paper 1 for a full description). 
We use these distributions to define time intervals of sim- 
ulated observations in our Monte Carlo analysis in order 
to determine our completeness in binary detections. 

This analysis shows that with our first two observa- 
tions we should detect 78% of the binaries with periods 
less than 10 3 days. If we include a third observation, 
this value increases significantly to 88%. If we include 
a fourth observation, our completeness only increases to 
90%. Thus to optimize our efficiency and completeness in 
detected binaries, we require three observations spanning 
roughly one year before classifying a star as either single 
or binary, confirming the results of[Mathicu (1983). We 
plot our binary detection completeness as functions of 
binary orbital elements in the solid lines in Figure [TJ and 
we use this analysis to correct our binary frequencies in 
Section [3] 

To estimate our incompleteness in binary orbital solu- 
tions, we simulate observations of the simulated binaries 
on the actual observation dates of each of the NGC 188 
binaries, respectively. For each combination of a simu- 
lated binary with a set of true observation dates, we at- 
tempt to derive an orbital solution using an automated 
version of the same code used in Paper 2. This code it- 
eratively minimizes the residuals between the RV data 
and the orbital fit to return the best fitting orbital so- 
lution. An acceptable orbital solution has RVs covering 
at least one period, errors on P, e and the orbital am- 
plitude K of less than 30% of the derived values, respec- 
tively, an RMS residual velocity of less than 1 km a~ l , 
and a range in RVs covering at least 75% of the orbital 
amplitude (mainly applicable for highly eccentric bina- 



ries). The resulting orbits define our completeness in 
orbital solutions, and each binary in NGC 188 produces 
a unique completeness curve in period, eccentricity, etc. 
The completeness resulting from all combinations of ac- 
tual observing dates and simulated binaries is plotted in 
the dashed lines in Figure [T] These histograms show the 
fraction of all binaries (within the specified period lim- 
its) for which we expect to derive orbital solutions, and 
hence obtain the orbital parameters. We use this incom- 
pleteness analysis to correct the observed distributions 
in period, eccentricity, secondary mass and mass ratio. 

Our completeness in detected binaries and orbital solu- 
tions decreases dramatically with increasing period (see 
the left panel of Figure [1} . Longer-period binaries gen- 
erally have a decreased amplitude of RV variations and 
require a longer time baseline of observations in order 
to detect and fully characterize the orbit. We detect 
88% of binaries with periods less than 1000 days, 78% of 
binaries with periods less than 3000 days, and 63% of bi- 
naries with periods less than 10 4 days. Only a negligible 
fraction of binaries with periods greater than 10 4 days is 
detectable with our current data. We therefore use the 
completeness in detected binaries with P < 10 4 days in 
our analysis of the binary frequency in Section |3l Only 
two of our MS binaries with orbital solutions have periods 
greater than 3000 days. Therefore for the following anal- 
ysis of the distributions of orbital parameters, we choose 
to limit our investigation to binaries with P < 3000 days. 

We show our completeness as a function of eccentric- 
ity fo r binaries with P r7 > r . < P < 3000 days (P C irc = 15.0 
days, iMathieu et al.l 120041 ) in the middle panel of Fig- 
ure [TJ The completeness in detected binaries is much 
less sensitive to orbital eccentricity than it is to period. 
However, our completeness in orbital solutions drops off 
rapidly at the highest eccentricities. High eccentricity bi- 
naries require well-placed observations to define the or- 
bital amplitude and eccentricity, greatly increasing the 
difficulty in deriving reliable orbital solutions. 

In the right panel of FigureQ] we plot our completeness 
in the mass ratio (q = MijM\) for binaries with P < 
3000 days. Because we use a single primary mass of 1 
Mq in this Monte Carlo analysis, the completeness in 
secondary mass is identical to the completeness in mass 
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ratio shown here. Therefore, this curve also shows our 
decreasing completeness with decreasing secondary mass, 
as expected from spectroscopic surveys of this type. 

The completeness curves in Figure [1] are not indepen- 
dent. Our detectability is most sensitive to period, since 
shorter-period binaries generally have higher-amplitude 
orbits. Hence this first bin in period shows the highest 
completeness, but the completeness is still less than 100% 
due to short-period binaries with very high eccentricity 
and/or very low- mass secondaries. Likewise no bin in ec- 
centricity or mass ratio shows 100% completeness due in 
part to long-period binaries found at all mass ratios and 
eccentricities. 

We use this analysis to correct our observed MS and 
giant binary frequencies and distributions of orbital pa- 
rameters for incompleteness. However we do not apply a 
completeness correction to the BS sample since we have 
generally observed the BSs more frequently than the ma- 
jority of the "normal" cluster stars, and the distributions 
of orbital parameters for the BSs may be significantly dif- 
ferent than the field dwarfs and subdwarfs surveyed by 
RIO. There is only one detected binary BS for which we 
do not have an orbital solution (WOCS ID 8104). 4 The 
orbital period for this binary is longer than the base- 
line of our observations. This system is accounted for in 
our analysis of the BS binary frequency, but not in the 
distributions of BS orbital parameters. 

Finally we note that there are only five non-RV- 
variable BSs in our sample. If we coarsely estimate their 
primary masses by t heir proximity to evolutionary tracks 
(jMarigo et al.1 12008) , and use these primary masses along 
with the true observations of these BSs in our Monte 
Carlo analysis (simplistically applying the same distribu- 
tions of orbital parameters as observed for field dwarfs), 
we find only a 24% chance of having an undetected bi- 
nary with a period < 10 4 days in this sample. Moreover 
we would only expect to have missed at most one BS 
binary in our data with P < 10 4 days. Therefore any 
incompleteness correction must necessarily be minimal. 

3. HARD-BINARY FREQUENCY 

A hard binary is defined as having an internal energy 
that is much greater than the energy of the relative mo- 
tion, with respect to the binary, of a sing le star movin g 
at the velocity disper sion of the c l uster (|Heggid fl975l ). 
Using the method of iGeller et al.1 (|2010| ) to correct for 
our measurement precision and undetected binaries, we 
derive a one-dimensional velocity dispersion in NGC 188 
of 0.41 ±g;5g km s" 1 for MS single members of NGC 188. 
(This value is consistent with that found in Paper 1.) 
Therefore, for solar-type stars in NGC 188, binaries with 
periods of less than ^10 6 days are hard binaries. As 
shown in Section [2] nearly all of our detectable binaries 
have periods less than 10 4 days. Therefore we detect 
only hard binaries, and in the following, we will use the 

4 There is also one proper-motion member in the BS region of the 
color-magnitude diagram for which we have been unable to derive 
a reliable RV membership due to rapid rotation and RV variability 
(WOCS ID 4230). The mean RV is outside of the cluster distri- 
bution, but we cannot be certain about membership without an 
orbital solution (and therefore a center-of-mass RV). We do not 
include this system in the following analysis due to it's uncertain 
membership; although we m ention it here because it is an X-ray 
source (S27, Gondoin 2005) and therefore warrants further obser- 
vations. 



term "hard binary" to refer specifically to the detectable 
hard binaries in NGC 188. These hard binaries are cat- 
alysts for the producti on of anomalous stars like BSs 
and X-ray sources (e.g., Bailvn 1995 1: iHurlev et alJl2005t 
llvanova et al.l 120061 ) , and as such the hard-binary fre- 
quency may directly influence the number of such anoma- 
lous stars produced within an open cluster. In this sec- 
tion, we compare the hard-binary frequencies of the NGC 
188 MS, giant and BS populations. 

We find 375 MS members including 90 MS binary 
members 5 . This observed MS hard-binary frequency is 
slightly inflated due to the faint magnitude limit defin- 
ing our stellar sample. The combined light of two MS 
stars in a binary system can be up to 0.75 magnitudes 
brighter than the primary star alone. Thus, the faint 
magnitude limit of our sample includes binary systems 
with primary stars that would lie below our magnitude 
limit. Here we discard the four binaries that, after cor- 
recting for secondary light, have primaries below the faint 
limit. We find an observed V magnitude corrected MS 
hard-binary frequency of 23 ± 2 % (86/371). As shown 
in Section [2] we detect 63% of the MS binaries with pe- 
riods of less than 10 4 days. The Monte Carlo analysis 
assumes that these binaries are detectable from the first 
three observations alone. However continued observa- 
tions of some stars in our sample have allowed us to de- 
tect certain (generally long-period) binaries after further 
observations, that did not show significant RV variability 
initially and therefore were not detected as binaries af- 
ter only the first three observations. Therefore, in order 
to correct this observed MS binary frequency for our in- 
completeness, we consider only the 67 binaries that show 
significant RV variability (i.e., e/i > 4, as in Paper 1) in 
their first three observations. When this sample is cor- 
rected for incompleteness, we find that the true MS bi- 
nary frequency among binaries with periods of less than 
1 4 days is 29 ± 3 %. 

iSollima et all d2010h estimate a MS core binary fre- 
quency for NGC 188 of 58.2 ± 13.5 % (over all mass ra- 
tios) through analysis of the CMD. To compare with this 
result, we first limit our sample to cover a similar spa- 
tial extent as that of ISollima et al.l (|2010l ). Specifically, 
we select all stars within 17 arcminute s from the clus- 
ter center. (We note that ISollima et al.1 ()2010l ) included 
lower-mass primaries than are included in our RV survey, 
but we do not attempt to correct for this difference in pri- 
mary mass range.) Within our limited sample, we find 
a solar-type MS binary frequency of 30 ± 4 % for bina- 
ries with periods < 10 4 days. If the NGC 188 solar-type 
MS binary period distribution is consistent with that of 
similar binaries in the Galactic field, then we can use the 
R10 field binary period distribution to extend this NGC 
188 binary frequency out to the hard-soft boundary of 
the cluster, of ~ 10 6 days. This analysis yields a MS 
binary frequency of 62 ± 8%, within the given spatial 
domain and out to the cluster hard-soft boundary, in ex- 
cellent agreement w ith the binary frequency derived by 
ISollima et all pTol ). 

The NGC 188 hard-binary frequency also agrees with 
values derived for other open clusters of a wide range of 

5 Here, as in Paper 1, we use our binary member (BM) and 
binary likely member (BLM) classes to define our binary member 
sample. 
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ages. iMathieu et"aH (|1990t ) find M67 (4 Gyr) to have a 
frequency of binaries with periods less than 10 3 days of 
9% to 15%. Matching this period range, we derive a bi- 
nary frequency of 22 ± 3 % out to a period of 10 3 days 
in NGC 188, somewhat higher than M67, but only at 
the margin of significance. The similar WOCS study 
of the intermediate-ag ed (2.5 Gyr) open cluster NGC 
6819 (Hole et al. 2009) report an observed binary fre- 
quency of 17%, all of which have periods less than 10 4 
days. We note that this observed binary frequency in 
NGC 6819 has not yet been corrected for incomplete- 
ness. The WOCS study of the young (150 Myr) open 
cluster M35 finds an incompleteness-corrected MS binary 
frequ ency of 24 ± 3% for binaries wi th periods < 10 4 
days (iGeller et alj|2010h . Additionally. iMermilliod et al.l 
(2008) hnd a spectroscopic binary frequency for stars of 
FGK spectral types of 20% in Blanco 1 (-100 Myr). 
Their observations have a similar time baseline to our 
NGC 188 survey though with fewer epochs, suggesting 
a slightly lower detection frequency as found here. Re- 
markably, these clusters span nearly 7 Gyr in age without 
a dramatic change in binary frequency. 

Out of the 454 primary stars in the R10 sample, there 
are 87 companions with periods less than 10 4 days, yield- 
ing a field binary frequency within this period range of 
19 ± 2 %. The NGC 188 hard-binary frequency out to 
the same period cutoff is marginally higher than that 
of the field. This higher binary frequency in NGC 188 
may be a dynamical signature, as TV-body models pre- 
dict that single stars will be preferentially lost from a 
cluster through evaporation a s compared to the gener- 
ally higher total mass binaries (|Hurlev et a l. 2005) . Thus 
NGC 188 may have started with a lower solar-type binary 
frequency, perhaps very similar to that of the field, that 
increased throughout the cluster lifetime. Furthermore, 
the hard-binary frequency in the halo of the cluster (see 
Figure [2]) , where the dynamical relaxation time is long 
and stellar densities are low, is in good agreement with 
that of the field. 

We find 70 giants, 19 of which are in binaries, resulting 
in an observed giant hard-binary frequency of 27 ± 6 %. 
(See Paper 1 for our definition of the giant region in the 
color-magnitude diagram, CMD.) 14 giant binaries would 
be detected after the first 3 observations alone, and are 
therefore applicable to our incompleteness correction. As 
a star evolves along the giant branch, its increasing phys- 
ical radius will preclude more and mo re companions at 
close separations and short periods. Mermillio d et al.l 
(2007) find a minimum period of 41.5 days in their sam- 
ple of red giant binaries in open clusters. Using this as 
our minimum period in the incompleteness analysis, we 
find a giant binary frequency of 37 ± 10 % for binaries 
with P < 10 4 days. The minimum period for the gi- 
ants in our NGC 188 sample is —11.5 days. Using this 
as the minimum period results in a binary frequency of 
34 ± 9 % for binaries with P < 10 4 days. The giant hard- 
binary frequency derived using either minimum period is 
consistent with that of the MS within the same period 
range, and with too large an error to see evolutionary 
differences in the frequency. 

In Figure [2] we plot the incompleteness corrected hard- 
binary frequency for the MS and giant populations as 
a function of projected radius from the cluster center. 
Again, we remove the four faint MS binaries with pri- 
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Fig. 2. — Incompleteness corrected hard-binary frequency as 
a function of projected radius, including the main-sequence and 
giant populations. We remove faint main-sequence binaries with 
primary stars that would likely fall below our observational limit 
of V = 16.5. The vertical error bars show the one sigma Poisson 
errors on the binary frequencies. The horizontal bars show the bin 
sizes of 5, 10 and 15 arcminutes, respectively. For main-sequence 
binaries we use our standard incompleteness correction, and for 
the giant binaries we apply a minimum period of 11.5 days in our 
incompleteness correction. We observe a central concentration of 
binaries as compared to the single stars in the cluster. 

maries that likely fall below our V magnitude limit. We 
exclude the BSs because their spatial distribution may 
be tied to their formation mechanism(s) rather than dy- 
namical (energy equipartition) effects. We divide our 
data into three bins of increasing size, covering 5, 10 and 
15 arcminutes, respectively, and denoted by the horizon- 
tal bars associated with each data point. As shown in 
Figure [5] the binary frequency clearly decreases with ra- 
dius from the cluster center. The binary frequency in the 
core is significantly higher than that of the halo; the first 
and last bins in Figure [5] are distinct with >99% con- 
fidence. This increase in binary frequency reflects the 
central concentration of the solar-type binaries with re- 
spect to the solar- type single stars in NGC 188 (as we 
also show in Paper 1). 

In marked contrast to the MS and giant populations, 
IMathieu k Gelled (|2009t ) find 16 of the 21 BSs to have bi- 
nary companions. (In Paper 1 we show the region on the 
CMD from which we have defined our BS population.) 
These 16 BSs are secure binaries as we have derived or- 
bital solutions for all but one; the remaining, 8104, is 
unquestionably a binary, with a longer period than our 
time baseline of observations for that star. Thus, we find 
a BS hard-binary frequency of 76 ± 19 %, more than 
three times and distinct with >99% c onfidence from the 
obser ved MS hard-binary frequency (|Mathieu k Gellerl 
2009). Furthermore, the five BSs that appear to be single 
may prove to be in longer-period binaries that currently 
lie outside of our spectroscopic dete c tion l imits. 

As discussed in IMathieu k Gelled (|2009f ). the large BS 
binary frequency is not sensitive to our selection criteria. 
For example, if we take only bright BSs (V < 15), the 
binary frequency is 76 ± 21 % (13/17). Alternatively if 
we include only those BSs with proper-motion member- 
ship probabilities >90% the binary frequency is 81 ±23% 
(13/16). 

We also note that there are two stars that lie above 
the MS turnoff and just fainter than our cutoff for BS 
selection (IDs 5101 and 8899). Our selection criteria ex- 
cludes any stars that could be members of normal MS 
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binaries. Thus, normal stars found just fainter than our 
BS region of the CMD are expected to be double-lined 
(SB2) binaries. To date we have not detected the flux 
from any companions in our spectra for either of these 
two stars, and their RVs do not vary above our thresh- 
old for binary detection. Therefore both of these stars 
are considered single cluster members (SM, see Paper 1), 
and perhaps should be included in the BS sample. If 
we include these two stars in the sample the BS binary 
frequency is 70 ± 17 % (16/23), still significantly higher 
than the observed MS binary frequency. 

In Paper 1 we find that the BSs appear to occupy a bi- 
modal spatial distribution, with a centrally concentrated 
population of 14 BSs and a halo population of 7 BSs. We 
compare the binary frequencies of the core and halo BS 
populations, divided at 10 arcmin from the cluster cen- 
ter, and find the core BSs to have a binary frequency of 
79 ± 24 % (11/14) and the halo group to have a binary 
frequency of 71 ± 32 % (5/7). The binary frequencies of 
these two subsets of NGC 188 BSs agree to within their 
large uncertainties. 

4. THE E - LOG(P) DIAGRAM 

Tidal circularization for the NGC 188 MS bi naries 
has been well studi ed by iMathieu et al.l (|2004l) and 
iMeibom fc Mathieul (|2005f ). Here we reintroduce the 
eccentricity versus log period (e — log(P)) diagram to 
present our new data, and to compare the MS, giant and 
BS distributions. We show the diagram in Figure [3] for 
the MS (top), giant (middle) and BS (bottom) popula- 
tions. Here, as for all such figures in this paper, we differ- 
entiate the MS, giant and BS binaries by plotting them 
in light gray, dark gray and black (green, red and blue 
in the online version), respectively. For reference, in the 
MS plot we show with the dotted line the circularization 
cutoff period, P C i rc = 15.0 days, found by M athieu et al.l 
(2004) and consist ent with the circularization period of 
14.5 days found bv IMeibom fc Mathieul (|2005| ). 

The MS and giant e — log(P) diagrams are of a sim- 
ilar form to those of solar-type binaries in the Galactic 
field (R10) and open clust ers of a wide range in age (see 
IMeibom fc Mathieul l2005t ) . The binaries are distributed 
across a large range in periods, with the short-period 
binaries (P < P C irc) having mainly circular orbits and 
longer-period binaries having a range in non-zero eccen- 
tricities. A two-dimension al Kolmogorov-Smirnov (K-S) 
test (using the technique of lFasano fc Fra nccschinii ll987l) 
comparing the MS and giant distributions reveals no 
statistically significant difference, although the expected 
lack of short-period orbits among the giants is present. 

However the BS binaries have a significantly different 
e — log(P) distribution than the MS (as also discussed in 
IMathieu fc GeHerll2009D . The majority of the BS binaries 
have periods on the order of 10 3 days. A two-dimensions 
K-S test results in 99% confidence that the MS and BS bi- 
naries are drawn from different parent populations. This 
concentration of BS binaries n ear periods of 10 3 days is 
also observed in the M67 BSs (|Lathamll2007f) as w ell as 
in metal-poor field BS binaries (|Carnev et al.ll2005l ) , and 
may provide significant insights into the origins of the 
BSs. 

The mean eccentricity for long-period (P > 500 day) 
BS binaries is 0.28 ± 0.06, while the mean eccentricity 
for similar MS binaries is significantly higher, at 0.44 
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Fig. 3. — Eccentricity plotted against the logarithm of the pe- 
riod for the main-sequence (top), giant (middle) and blue straggler 
(bottom) populations. For the main-sequence binar ies, we plot the 
tidal circularization cutoff period P c i rc = 15.0 days l IMathieu et al.l 
2004) with the dotted line. For the blue straggler binaries, we iden- 
tify the core (triangles) and halo (squares) populations. 



± 0.06 (jMathieu fc Geller! 12009] ) . Furthermore, three of 
these twelve long-period BS binaries have eccentricities 
within two standard deviations of zero, while none of the 
long-period MS binaries have such low eccentricities. 

Additionally, in the bottom panel of Figure [3] we di- 
vide our BS binaries to compare the core and halo BS 
populations (again divided at 10 arcmin from the cluster 
center). We plot the eleven core BS binaries with trian- 
gles and the four halo BS binaries with squares. We find 
no statistically significant distinction in either one- or 
two-dimensional K-S tests comparing the period and ec- 
centricity distributions of these two populations, though 
these results are drawn from small sample sizes. 

Finally, we note that none of our binaries with orbital 
solutions have periods of less than one day. Indeed our 
shortest-period MS binary has a period of 3.8 days, and 
our shortest-period giant binary has a period of 11.5 
days. We can understand this lack of short-period gi- 
ant binaries as an evolutionary effect due to increasing 
radius. For two solar-mass stars, a period of 10 days 
corresponds to an orbital semi-major axis of ^25 R©, 
which is smaller than the radii of the N GC 188 giants 
(jvan Belle et alJll999t iGirardi et al.ll2002t ). 

On the other hand the R10 period distribution 
predicts that we should observe ~6-8 MS binaries 
with periods below 4 days. It is well known that 
NGC 188 contains a large numbe r of W UMa sys- 
tems (e.g. iBaliunas fc Guinanl Il985l; IZhang et"ail 120021: 
iKafka k Honevcuttj 12003( 1 There are four known W 



WOCS. The Hard-Binary Population of NGC 188 



7 



UMa's within our sample that have non-zero proper- 
motion memberships, each having a photometric period 
of less than one day. We have been unable to obtain kine- 
matic orbital solutions for these systems, as their rapid 
rotation hinders our ability to derive precise RVs. This 
lack of detached short-period binaries combined with the 
high frequency of W UMa systems in NGC 188 suggests 
that we may be observing an evolutionary depletion of 
close binaries. Perhaps this is a result of orbital angular 
momentum loss from magnetic braking in stellar winds 
that has slow ly decreased the semi-majo r axes of the clos- 
est binaries ((Baliunas fc Guinanl 1 19851 ) . Alternatively, 
short-period binaries in triple systems may have been 
driven t ogether through the Kozai mechanism and tidal 
friction (jEggleton fc Kisseleva-Eggletonl2006l) or dynam- 
ical encounters. (We have not yet detected tertiary com- 
panions to any of the W UMa's in our sample, though 
kinematic detection of tertiaries is impeded by the rapid 
rotation of the W UMa.) 

5. ECCENTRICITY DISTRIBUTION 

Here we discuss the eccentricity distribution for bina- 
ries with periods between P C j rc and 3 000 days, where 
P circ = 15.0 days (jMathieu et al.ll2004h . In Figure H we 
plot the eccentricity distributions of NGC 188 binaries 
in both histogram and cumulative distribution form. We 
plot the detected binaries with orbital solutions in filled 
histograms. The incompleteness-corrected distributions 
(for MS and giant binaries) are plotted with the solid- 
lined histograms above the observed distributions. We 
choose to bin the MS data with bins of 0.1 and the gi- 
ant and BS data with bins of 0.15 to reduce the scatter 
and larger errors associated with the relatively smaller 
number of giant and BS binaries. In the bottom panel 
of Figure HI we plot the cumulative distributions of these 
NGC 188 binary populations. 

We correct the MS and giant distributions for incom- 
pleteness by applying the results of Section [2] (Again, 
we do not account for evolutionary effects in our Monte 
Carlo incompleteness study; e.g., to account for the in- 
creased radius of a giant star). As discussed in Section [21 
we do not attempt to correct the BS distributions for in- 
completeness as the masses and distributions of orbital 
parameters for these primaries are likely significantly dif- 
ferent from the 1 M Q primaries and field dwarf binary 
orbital parameters used in our incompleteness analysis. 
Again note that there are only six additional BSs in NGC 
188, only one of which is detected as a binary. Thus any 
incompleteness correction for the BS binaries must be 
minimal. 

We also show a Gaussian fit to the R10 field eccentric- 
ity distribution and a thermal distribution (/(e) = 2e) 
as the dashed and dotted lines, respectively, both nor- 
malized by the respective sample sizes using the R10 
field binary frequency of 11.6% in this period range. In 
the cumulative distribution plot, we show the R10 and 
thermal distributions modified by our observational com- 
pleteness (in the corresponding line styles). To produce 
these curves we multiply the R10 and thermal distribu- 
tions, respectively, by our incompleteness distributions 
in orbital solutions (from Figure [lj , effectively limiting 
these distributions to reflect our completeness level as 
a function of eccentricity. These completeness-corrected 
distributions are then converted into the cumulative dis- 
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Fig. 4. — Eccentricity distributions for NGC 188 binaries with 
Pcirc < P < 3000 days. We show the observed main-sequence, 
giant and blue straggler binaries in light gray, dark gray and black 
(green, red and blue in the online version), respectively. The top 
three plots show the distributions of eccentricity in histogram form. 
We correct the main-sequence and giant distributions for our in- 
completeness and plot the corrected histograms with the solid lines 
above the observed distribu tions. In each hi s togra m plot, we also 
show a Gaussian fit to the Raghavan ct al. (2010) field distribu- 
tion in the dashed line and a thermal distribution in the dotted 
line, both normalized to the respective sample size using the field 
binary frequency within this period range of 11.6%. The bottom 
plot shows the cumulative distribution functions for the three pop- 
ulations in the same gray-scale (color in the online version) coding 
as for the histograms. For clarity, we plot the main-sequence bina- 
ries with circles, giant binaries with diamonds, and we distinguish 
the core and halo blue straggler populations with the triangles and 
squares, respectively. As a comparison, we also plot the field and 
thermal distributions corrected for our observational completeness 
in orbital solutions in the dashed and dotted lines, respectively. 

tributions shown in Figure [4] 

The NGC 188 MS eccentricity distribution has a 
roughly Gaussian form centered on e — 0.35. We find no 
distinction between the MS NGC 188 eccentricity distri- 
bution and either the Gaussian fit to the R10 eccentricity 
distribution or a flat eccentricity distribution. The NGC 
188 MS eccentricity distribution also agrees with that of 
similar samples in open clusters of a wide range in age. 
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Fig. 5. — Eccentricity as a function of projected radius from the 
cluster center. We plot the main-sequence binaries in light-gray 
circles, giant binaries in dark-gray diamonds and blue straggler bi- 
naries in black triangles (core) and squares (halo). (We use green, 
red and blue for the main-sequence, giant and blue straggler bi- 
naries, respectively, in the online version.) Note that we have not 
derived orbits for any binaries with e > 0.5 outside of 5 core radii 
from the cluster center. When we divide the NGC 188 binaries at 
5 core radii (~ 13 arcminutes), we find the inner binaries to have a 
distribution shifted to higher eccentricities at the 97% confidence 
level. 

The mean eccentricity of similar binaries in both M67 (4 
Gyr) and the Pleiades (150 M yr) is -0.4 (jMathieu et al.l 
119901 : iMermilliod et all 11993 ). Likewise, the eccentric- 
ity distribution of M35 (—150 Myr) is well approximated 
by a Gaussian with a mean eccentricity of — 0.3 - 0.4 
(jMathieul 120081 ). Conversely a K-S test shows that the 
NGC 188 MS eccentricity distribution is not drawn from 
a thermal distribution at the >99% confidence level. 

We find no statistically significant distinction between 
the giant and MS eccentricity distributions. However a 
K-S test comparing the BS and MS eccentricity shows 
the two distributions to be distinct at the 97% confi- 
dence level, with the BS distribution shifted to lower 
eccentricities. This res ult agrees with that found by 
iMathieu fc Gelled (|2009f) and also discussed in Sectional 
where we show that the BS binaries with P > 500 days 
have significantly lower eccentricities than MS binaries 
of similar periods. 

We also investigate the dependence of eccentricity on 
projected radius from the cluster center. In Figure [5] we 
show all binaries with orbital solutions, no longer lim- 
iting by period. As is clear from this figure, we have 
not derived orbital solutions for any high eccentricity 
(e > 0.5) binaries at r > 5 core radii (r > 13 arcmin- 
utes), while we do have such binaries in the inner-cluster 
region. If we divide our sample of all binaries at 5 core 
radii, we find the inner binaries to have a mean eccen- 
tricity of 0.31 ± 0.03, and the outer binaries to have a 
mean eccentricity of 0.16 ± 0.04. These two mean val- 
ues can be distinguished at the 97% confidence level by 
a Student's T-test. Likewise a K-S test shows that the 
inner binaries have a distribution shifted to higher ec- 
centricities at the 97% confidence level. This distinction 
may be the result of dynamical encounters that are ex- 
pected to occur most frequency in the denser core of a 
cluster. In star clusters, fly- by interactions are the most 
common type of encounter (jBacon et al.lll996l ). and can 
per turb a binary causing an increase in orbital eccentric- 
ity (|Heggie fc R asio 1996). Furthermore close three- and 
four-body interactions involving lower eccentricity bina- 
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Fig. 6. — Period distributions for the NGC 188 binaries with 
P < 3000 days. We show the observed main-sequence, giant and 
blue straggler binaries in light gray, dark gray and black (green, red 
and blue in the online version), respectively. The top three plots 
show the period distributions in histogram form, while the bot- 
tom plot shows the cumulative distributions. The last bin extends 
to periods of 3000 days and is normalized to reflect the different 
bin size. We correct the main-sequence and giant binary period 
distributions for incompleteness and plot the corrected histograms 
with the solid line, above th e observed distribution s. In each his- 
togram plot, we also show the Raghavan et al. (2010) Galactic field 
log-normal distribution (dashed line) normalized to the respective 
sample sizes using the field binary frequency within this period 
range of 14.1%. In the bottom plot, we show the cumulative distri- 
butions of our three NGC 188 samples in the same gray-scale (color 
in the online version) coding as for the histograms. Again, we plot 
the main-sequence binaries in circles, giant binaries in diamonds 
and blue stragglers binaries in triangles (core) and squares (halo). 
Here we plot the field distribution corrected for our observational 
completeness in orbital solutions with the dashed line. 



ries have a high likelihood of producing high-eccentricit y 
products (|Fregeau et al J 120041: iGiersz fc Spurzemll2003f) . 
However, we caution the reader that our completeness 
in orbital solutions drops off dramatically for very high 
eccentricity binaries (Figure [TJ . 
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Fig. 7. — The logarithm of the period as a function of projected 
radius from the cluster center. We plot the main-sequence bina- 
ries in light-gray circles, giant binaries in dark-gray diamonds and 
blue straggler binaries in black triangles (core) and squares (halo). 
(We use green, red and blue for the main-sequence, giant and blue 
straggler binaries, respectively, in the online version.) We find no 
trend in period with radius. (Note the errors in period, when plot- 
ted in the log, lie within the symbols used for plotting the data 
points.) 

6. PERIOD DISTRIBUTION 

In Figure |B] we plot our observed MS, giant and BS 
period distributions in both histogram and cumulative 
distribution form for binaries with P < 3000 days. Addi- 
tionally we plot the R10 period distribution (dashed line) 
normalized to the respective sample sizes using the R10 
binary frequency within this period range of 14.1%. For 
the MS and giant distributions we use our incompleteness 
study of Section [5] to correct our observed values. In the 
bottom panel of Figure [6l we plot the cumulative distri- 
butions of the NGC 188 binaries as compared to the R10 
distribution modified by our observational completeness 
in the same manner as in Figure [4] 

The incompleteness-corrected NGC 188 MS log-period 
distribution rises towards our completeness limit, in good 
agreement with the log-normal distribution found by R10 
for solar-type field binaries. A K-S test comparing the 
MS NGC 188 and R10 cumulative distributions results 
in a 67% probability that they are drawn from different 
parent distributions. This value is as large as 67% due to 
the lack of observed binaries with periods less than ^4 
days, as discussed in Section 21 Thus we find that the 
MS binary period distribution is indistinguishable from 
the R10 field distribution over our period range. 

A similar analysis returns >99% confidence that the 
present day NGC 188 MS binary population is not drawn 
from a flat distribution in log(P), a common choice for 
initializing binaries in iV-body models. Additionally the 
flat period distribution is not supported by obs ervations 
of th e young (~150 Myr) open cluster M35 (|Mathieul 
l2008h . 

A K-S test comparing the cumulative distributions of 
the MS and BS binary populations results in >99% con- 
fidence that these t wo samples are drawn fr om distinct 
parent populations (jMathieu fc Gellerll2009l ). The ma- 
jority of the BS binaries are found at periods near 1000 
days, as is also clear from Figure [31 

A similar test comparing the giant and MS binaries 
results in an 88% probability that they are drawn from 
different parent populations. This marginal distinction 
between the giants and MS binaries is partly due to the 



lack of short-period giant binaries. This loss of short- 
period binaries is expected as the primary evolves off the 
MS and expands to become a red giant. If we exclude 
MS binaries with P < 11.5 days (the minimum orbital 
period found amongst the giant binaries in NGC 188), 
the probability that the MS and giant distributions are 
different is reduced to <50%. 

We find no dependence in period with projected radius 
from the cluster center (Figure [7]). The core and halo BS 
binary populations are seen clearly in Figure [7] There is 
no distinction between the period distributions of these 
two subsets of the BSs, even to the extent of each having 
one short-period SB2 BS binary. 

7. SOLAR- TYPE MAIN-SEQUENCE SECONDARY-MASS 
AND MASS-RATIO DISTRIBUTIONS FOR HARD 
BINARIES 

The secondary-mass and mass-ratio distributions for 
most binary populations are less well known than, for in- 
stance, the period and eccentricity distributions. Kine- 
matic orbital solutions alone do not reveal the compo- 
nent masses for spectroscopic binaries. For SB2 binaries 
we derive the mass ratio directly; however the individ- 
ual masses obtained from the kinematic orbital solutions 
are dependent on the inclination angle. For single-lined 
(SB1) binaries, we only derive the mass function in which 
the primary and secondary masses are entangled and also 
dependent on the inclination. 

In Paper 2 we used a photometric deconvolution tech- 
nique to estimate primary masses for all of the MS bi- 
naries with orbital solutions. Here we use these pri- 
mary mass estimates and the dynamical mass functions 
within a statistical algorithm (Mazeh & Goldbcrj |1992l) 
to derive the secondary-mass and mass-ratio distribu- 
tions. In Figure [5] we plot the resulting MS mass-ratio 
and secondary-mass distributions for the NGC 188 bi- 
naries. As in previous figures, all detected binaries with 
orbital solutions are shown in the filled histograms, while 
the incompleteness corrected distributions arc plotted in 
the thick solid lines above the observed distributions. 
The horizontally-hatched regions in Figure [8] show the 
SB2 binaries and the vertically-hatched regions show the 
SB1 binaries. 

For the SB1 bina r ies we use the iterative algorithm of 
iMazeh ^G oldberg (|1992h to convert our distributions of 
mass functions into distributions of mass ratio and sec- 
ondary mass, using our primary mass estimates and as- 
suming the inclination angles to be randomly distributed. 
Note that this technique (as well as other statistical in- 
version methods) does not reveal the secondary mass or 
mass ratio for a given binary, but instead derives the 
respective distributions for a sample of binaries. 

We correct the SB1 binaries for our incompleteness us- 
ing the analysis of Section [2] after applying the statistical 
algorithm to derive the observed distributions. We do 
not use this same incompleteness correction for the SB 2 
binaries as these systems (once detected) are given the 
highest priority in our observations. As such, we have 
only one confirmed SB2 binary without an orbital so- 
lution (ID 9119), and we add in this binary by hand. 
The turnoff in NGC 188 is at ~1.1 M©, and our incom- 
pleteness correction covers binaries with P < 3000 days 
that have stellar companions down to 0.08 M (exclud- 
ing brown dwarfs and other sub-stellar objects). 
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Fig. 8. — Secondary-mass (left) and mass-ratio (right) distributions for the NGC 188 solar-type main-sequence hard binaries. The 
detected binaries with orbital solutions are plotted in the filled histograms, with SB2s in the horizontally- hatched histograms and S Bls 
in the vertically-hatched histograms. For the SB1 binaries we utilize the statistical inversion technique of Mazch & Goldberg (1992) to 
derive the secondary- mass and mass-ratio distributions from the observed mass functions. The incompleteness-corrected distributions are 
plotted as the thick solid lines above the observed distributions. The error bars are derived through a Monte Carlo analysis utilizing the 
Mazch & Goldberg (1992) algorithm to convert the Poisson errors on the mass function distribution into uncertainties on the secondary- 
mass and m ass-ratio distributions, respectively, and we show the 95% confidence intervals. For comparison, in the secondary-mass panel 
we show the Kroupa (2001) IMF, and in the mass-ratio panel we plot the distribution obtained from choosing random partners from this 
IMF, both with the solid curves. The curves are normalized to contain the same number of binaries as in our incompleteness-corrected 
NGC 188 sample. Finally we plot the secondary-mass distribution in the more familiar log-log format in the upper-right inset. We show 
this figure for illustrative purposes, but focus our statistical analysis on the observed mass-function distribution shown in Figure [9] 
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Fig. 9. — Cumulative mass-function distribution for the NGC 
188 solar-type main-sequence binaries. We show the observed mass 
functions with the black points, and circle the double-lined binaries 
(for reference). With the solid line we plot the ma ss function dis - 
tribution derived by drawing companions from the Kroupa (2001) 
IMF. To derive this curve we assume the inclination angles are dis- 
tributed isotropically and draw primary masses from a Gaussian 
fit to the observed NGC 188 main-sequence primary-mass distri- 
bution. This theoretical distribution is then corrected for our in- 
completeness in orbital solutions (out to P = 3000 days) in the 
same manner as Figures [4] and [6] A K-S test shows that the ob- 
served mass-function distribution is statistically indistinguishable 
from the IMF. 



In the secon dary-mas s plot , we also show the "univer- 
sal IMF" from iKroupal (|2001h with the solid curve. In 
the mass-ratio plot, we show the distribution obtained 
from choosing random partners from the IMF in the solid 
curve, limiting the primary to be within our observable 
mass range and with the mass ratio always less than one. 
The curves are normalized to contain the same number 
of binaries as in our incompleteness corrected NGC 188 
sample^ 

The IKroupal ()2001l ) IMF was not intentionally derived 
to model secondary stars in binaries. Still the shape of 
the MS secondary-mass distribution is consistent with 



the IMF in that it rises towards lower masses. Like- 
wise the mass-ratio distribution shows a similar form to 
the curve derived by choosing random partners from this 
IMF. 

We also observe an excess of binaries with masses be- 
tween 0.85 and 1.0 M (and mass ratios > 0.9) as com- 
pared to the IMF. This excess of near equal-mass bina- 
ries remains if we instead use the mass functions for all 
binaries (including the SB2s) and derive the secondary- 
mass and mass-ra tio distributions for the full binary pop- 
ulation with the iMazeh fc Goldberg! (|1992fl algorithm; 
therefore this result is not sensitive to the method used 
here. This may be analogous to the abundance of "twins" 
obse rved in some studies, generally at shorter periods 
(e.g. lTokovinh]|2?j00t iFisher et al.ll2005l R10). We have 
investigated the mass-ratio distribution as a function of 
period, but find that the decreased sample sizes make 
these results highly uncertain. 

To check formally for consistency between the observed 
companion masses and the IMF we turn to the observed 
mass functions. In Figure^ we compare the observed cu- 
mulative mass-function distribution (black points, with 
SB2s circled) to a theoretical distribution (solid line) de- 
rived by drawing companions from the IMF, primary 
masses from a Gaussian fit to our observed NGC 188 MS 
primary mass distribution, and assuming an isotropic in- 
clination distribution. We then correct this theoretical 
distribution for our incompleteness in orbital solutions 
using the same method as in Figures |4] and [6] A K-S test 
shows that we cannot formally distinguish the observed 
mass-function distribution from the IMF. (We also note 
that the mass-function distribution cannot formally be 
distinguished from a distribution derived by choosing bi- 
naries from a uniform mass-ratio distribution.) 

There is some disagreement in the literature about 
the secondary-mass and mass-ratio distributions for field 
dwarfs. R10 find their sample of field dwarfs to have a 
uniform mass-ratio distribution between 0.2 < q < 0.95, 
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Fig. 10. — Secondary-mass (left) and mass-ratio (right) distributions for the NGC 188 blue straggler binaries. The SB2s are plotted with 
the horizontally-hatched histograms and the SBls are plotted wi th the vertically-hatc hed histograms. For SBls, the primary masses are 
estimated in relation to evolutionary tracks of Marigo ct al. (2008), and we then use the Mazeh & Goldberg (1992) algorithm to derive the 
secondary-mass and mass-ratio di stributions from the obse rved mass functions. The masses for the two SB2 blue stragglers are discussed 
in detail in Section[8] and see also Mathicu & Gcllcr (2009). The uncertainties shown here are derived in the same manner as in Figure 151 



with evidence for an abundance of twins at q > 0.95. 
iMazeh et al.l (|2003l) find a mass-ratio distribution that is 
formally consistent with a uniform distribution. They 
also note a possible rise in the distribution towards 
lower mass-ratio values, bu t wi th large uncert a inties . 
iDuauennov fe Mayor! (|199lD and IGoldberg et~aTl (|2003f ) 
find distribution s that rise towards lo wer-mass secon- 
daries. (Strictly IGoldberg et~all (|2003| ) find this result 
for disk binaries, and not for the halo binaries.) We con- 
clude that the companions to the NGC 188 MS binaries 
are consistent with those of the field dwarf binaries, but 
given the uncertainty in the literature about the field 
population this conclusion lacks precision. 

8. BLUE STRAGGLER SECONDARY-MASS AND 
MASS-RATIO DISTRIBUTIONS 

The masses and evolutionary states of the companion 
stars to the BSs in binaries are a powerful tool to investi- 
gate the origins of these systems. Knowledge of whether 
the secondaries are normal MS stars or white dwarfs 
(WDs) can be particularly useful in distinguishing be- 
tween possible formation mechanisms (|Geller fe Mathieul 
1 2 1 ll ) . For instance, the Case B or Case C mass transfer 
mec hanisms predic t a BS binary w ith a WD compan- 
ion (|McCreal TT964; Pac zvriskl 1 19711 ), while a merger of 
an in ner binary of a triple system ([Perets fe Fabrvckvl 
[III or a stellar collisi on that leaves a bound compan- 
ion ^Hurley et al.ll2005[ ) would more likely produce MS 
secondaries. 

In order to derive the observed BS secondary-mass and 
mass-ratio distributions, we must first estimate the BS 
primary masses. To do so we compare the observed mag- 
nitudes and col ors of the BSs to a coarse grid of evolu- 
tionary tracks (jMarigo et al.l 120081 ) extending from the 
ZAMS to 7 Gyr, and for stars with masses ranging from 
the cluster turnoff mass of 1.1 M© to twice this mass (in 
steps of 0.1 M Q ). For the SB1 BSs, the secondaries are at 
least 2.5 magnitudes (in V) fainter than the BS primaries 
(Paper 1). For a given BS, deconvolving a secondary 
star at this limiting magnitude from the combined light 
would result in a primary of < 0.1 V magnitudes fainter 
than the combined V magnitude. This minor difference 
in V magnitude is below the mass resolution of the evo- 



lutionary tracks (e.g., see Figure ITTT) . Therefore we do 
not attempt to deconvolve the secondary light from these 
BSs, and estimate their masses based on the proximity 
of the combined magnitude and color to the evolutionary 
tracks. We find the NGC 188 SB1 BSs to have masses 
between 1.2 - 1.6 Mq. 

We then use the Ma zeh fc Goldberg! (jl992) algorithm 
in the same manner as for the MS binaries to derive the 
BS binary secondary-mass and mass-ratio distributions 
shown in Figure [10] (with SBls in the vertically cross- 
hatched histograms and SB2s in the horizontally cross- 
hatched histograms). The BS SB1 secondary-mass distri- 
bution is narrow and peaks near a mass of ~0.5 Mq. This 
result, when considered in light of the periods near 1000 
days and low eccentricities, suggests an origin through 
Case C mass transfer, which predicts carbon-oxygen WD 
companions of masses between ~0.5 - 0.6 Mp, with or- 
bital periods of order 100 davs dMcC rca 1964; Pac zvfiskl 
[19711: iHurlev et al.l 120021: iChen fc Hanll2008D . However 
the secondary-mass distribution alone does not rule out 
MS companions, as we would not detect the flux from a 
0.5 M Q MS (or WD) star in our BS spectra due to the 
large di fference in luminosi t ies. W e discuss this result in 
detail in iGeller fc Mathieul (|201 lh and build upon these 
results in Section [H 

For the two short-period SB2 BS binaries, the sec- 
ondary adds significant light to the system, and thus for 
each SB2 BS we deconvolve the secondary from the com- 
bined light before estimating the primary mass. Addi- 
tionally we utilize a grid of synthetic spectra to estimate 
the effective temperatures of both components of these 
systems, which helps to further constrain the masses. We 
discuss these results briefly in iMathieu fc Gellerl ([2009D 
and expand upon them here. 

BS ID 7782 has a kinematic mass ratio of q = 1.005 ± 
0.013 (Paper 2). A spectral analysis yields a luminosity 
ratio 6 of 0.871 ± 0.024, and effective temperatures for 
the primary and secondary stars of 6500 K and 6325 K, 
respectively. Both stars are hotter than the MS turnoff 

6 Due to an update in our TODCOR analysis we fin d a slightly 
different luminosity ratio from Mathicu & Gcller (2009), who gave 
a value of 0.75 ± 0.01. 
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Fig. 11. — Color-magnitude diagram showing the loci of possible 
locations for the components of the blue straggler binary ID 5078. 
Cluster members are shown in the small black points, and 5078 is 
circled. We also show a ZAMS and a 7 Gyr isochrone in the thin 
solid lines, and evolutionary tracks for 1.2 Mq, 1.3 Mq, 1.4 Mq an d 
1.5 Mq stars in the dotted lines, respectively (Marigo ct al. 2008). 
The thick gray line shows the locus of potential secondary stars, 
and the thick black line shows the locus of potential primary 
stars. The derived effective temperature of 5850 K for the sec- 
ondary implies a mass of 1.02 Mq; we plot the location of this 
secondary star, and the location of the associated primary, with 
the light-gray crosses. However, given the kinematic mass ratio 
of q = 0.678 ± 0.009 and this secondary mass, the blue straggler 
mass would be 1.5 Mq. Yet the blue straggler is significantly less 
luminous than a normal 1.5 Mq star at any point in its evolution. 

effective temperature of 5900 K and more luminous than 
turnoff stars. Evidently 7782 comprises two BSs. The 
location of the system in the cluster halo suggests that a 
close dynamical encounter was involved in the formation 
of the binary, perhaps the final of several encounters near 
the cluster core that exchanged the BSs into the system 
or formed one of them collisionally and ejected the binary 
into the cluster halo. 

The differing temperatures and luminosities are unex- 
pected given the nearly identical masses. Possibly the 
two BSs did not form at the same time and thus have 
different evolutionary ages. We are able to model the 
system wit h two 1.25 Mp star s having an age difference 
of ^2 Gyr (jMarigo et al.ll2008h . although this solution is 
not unique given the present data. More detailed spec- 
tral analysis of the two stars is merited. 

Nonetheless, given the large bin sizes used in construct- 
ing the secondary-mass distribution (Figure [10]), we can 
confidently place the secondary of 7782 in the 1.2 - 1.4 
M Q bin. 

BS ID 5078 has a kinematic mass ratio of q = 0.678 ± 
0.009 (Paper 2). From our analysis of the spectra, we 
derive an effective temperature for the primary of 6500 
K and for the secondary of 5850 K. The effective tem- 
perature of the secondary is consistent with a star on 
the upper MS. Assuming the secondary is a MS star, 
we deconvolve the combined light using secondary light 



contributions derived from a 7 Gyr Mar igo et al.l (2008) 
isochrone. We limit the possible secondaries such that 
the effective temperature is within 250 K (~ 2a) of 5850 
K, and we limit the possible BS primaries such that the 
BS is on or redder than the ZAMS. The resulting loci of 
possible CMD locations for the primary and secondary 
are shown in the thick black and gray lines in Figure [Til 
respectively. 

We find a range in possible masses for the secondary 
of 5078 of 0.94 M to 1.08 M Q . The nominal effective 
temperature for the secondary of 5850 K corresponds to 
a 1.02 M Q MS star, and we use this secondary mass in 
Figure [10] The location of the associated BS primary 
lies on or near a 1.2 Mq evolutionary track. We note 
though that the effective temperature of 6500 K derived 
for the BS primary suggests a (B — V) color that is ~0.1 
mag bluer than inferred from this analysis. 

The range in secondary masses and the kinematic mass 
ratio yields a range for the BS primary star mass between 
1.39 M Q and 1.59 Mq. However, the locus of potential 
CMD locations for the primary clearly falls below such 
evolutionary tracks, and even the combined light of 5078 
is too faint. The standard stellar evolution tracks un- 
derestimate the dynamically informed mass of the BS 
primary star by 15% to 30%. 

This potential systematic error does not meaningfully 
effect the results shown in Figure [101 If we simply assign 
the same mass to each BS within the observed range of 
1.2 - 1.6 M , the resulting secondary-mass and mass- 
ratio distributions stay within the uncertainties on the 
distributions shown in Figure [TTj] Assigning all BSs a 
mass of 1.3 Mq, which is approximately the average BS 
mass that we derive, results in distributions that are 
nearly identical to those shown here. If we take the ex- 
treme case of setting all BSs to twice the turnoff mass 
(2.2 Mq), the general forms of the distributions remain 
the same, but the peak in the secondary-mass distribu- 
tion shifts to the 0.6 - 0.8 M Q bin. 

Again we do not attempt to correct these BS distribu- 
tions for our incompleteness and show only the BS bina- 
ries with orbital solutions. We remind the reader that 
there are only six additional NGC 188 BSs, only one of 
which is detected as a binary but does not yet have an 
orbital solution. 

9. COMPARISON TO A SOPHISTICATED W-BODY OPEN 
CLUSTER SIMULATION 

9.1. The Simulation 

The open cluster TV-b ody simulation of M67 by 
iHurlev et ail (|2001l |2005|) was able to track the evolu- 
tion and dynamics of realistic initial numbers of single 
stars and binaries, and included detailed prescriptions 
for mass loss, tidal processes, mass transfer, mergers and 
collisions. Through their model, Hurley et al. were able 
to match the observed M67 color-magnitude diagram re- 
markably well. They produced similar populations of 
anomalous stars as are observed in M67, including BSs. 
Importantly, they were able to track the dynamical evo- 
lution of these anomalous stars and pinpoint their origins 
with great detail, many of which involved interactions 
with binary stars. Simulations such as these have the 
power to revolutionize our understanding of the evolu- 
tion of a binary population and the origins of anomalous 
stars, like BSs, in star clusters. 
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Though Hurley et al. designed their simulations to 
model M67, their results can also be compared to simi- 
lar open clusters, such as NGC 188. Both NGC 188 and 
M67 are old open clusters at 4 and 7 Gyr respectively. 
The curre nt dynamical relaxa tion time for NGC 188 is 
180 Mvr (IChumak et al1l2010h and for M67 is 100 Myr 
(|Mathieu fc Lathamlll986l) . As shown in Section El NGC 
188 and M67 have comparable binary frequencies, and 
their binary orbital pe riod and eccen tricity distributions 
are of a similar form ([Latham! 12007). Furthermore, the 
BS popul ations of these two clu sters show strong simi- 
larities. lLatham fc Milond (|1996| ) observed ~60% of the 
M67 BS stars to be members of hard-binary systems, 
consistent with the 76 ± 19 % binary frequency amongst 
the NGC 188 BSs observed here. Also, the M67 BS bi- 
naries sho w a similar tendency towa rds periods of order 
1000 days (jLatham fc Miloneill996T ). Thus, in this sec- 
tion we comp are the NGC 188 h ard-binary population 
to that of the iHurlev et al.l ()2005l ) simulation. 

The simulation started with a 50% binary frequency 
over all periods and masses, in which binary orbital el- 
ements were chosen from a flat period distribution and 
a thermal eccentricity distribution. The pr imary stars 
were chosen to follow a iKroupa et al.l ()1991l ) IMF, with 
masses > 0.1 M , and the mass ratios for binaries were 
chosen from a uniform distribution in q. The bi naries 
were then modified according to the IKrou pa (1995) tidal 
evolution prescription, and the resulting initial stellar 
population was evolved for 4 Gyr. 

We compare the MS and giant populations of the sim- 
ulation at this final age to NGC 188. We only analyze 
MS and giant stars that would be within our observable 
period, spatial and magnitude domain. Specifically we 
consider only those binaries with P < 10 4 days to be 
detectable, and assume that we would only obtain or- 
bital solutions for binaries with P < 3000 days. Our 
analysis only includes the inner 17 pc (radius in projec- 
tion along the z axis) of the cluster. We use the NGC 
188 distance modulus o f (M — w)y=11.44 an d the red- 
dening E(B-V)=0.09 (jSaraiedini et al.H l9991 to define 
the lower-mass (faint) limit of the MS. The giants are all 
within our observable magnitude domain; thus we only 
limit the giants by radius from the cluster center. 

We choose not to impose magnitude or spatial limits 
on the BSs created within the simulation. There are 
only four BSs created within the simulation throughout 
the entire cluster lifetime that at some moment reside 
outside of our observable spatial region. Each of these 
were formed near the core, and later ejected from the 
cluster, spending only one time step (^60 Myr) outside 
of 17 pc (in projection) from the cluster center before 
traveling beyond the cluster's tidal radius. 

In the following analysis, we use an "integrated sam- 
ple" of all BSs that existed between 2 and 4 Gyr. The 
"integrated sample" comprises the union of every BS at 
each time step in the simulation during this range in time. 
There are 57 individual BSs present at some point during 
this interval, while the integrated sample contains multi- 
ple realizations of these BSs over these 33 time steps. In 
so doing, the integrated sample effectively weights results 
by the lifetime of each individual BS in a certain binary 
configuration. For example, there were 26 BSs created 
in binaries between 2-4 Gyr, 11 of which underwent an 
encounter after formation that significantly changed the 




1 2 

log (Period [days]) 



Fig . 12. — Period distribution for binaries in the IHurlev et af] 
(2005) simulation. We compare the main-sequence (top), giant 
(middle) and blue straggler (bottom) populations from the simu- 
lation, using the 4 Gyr main-sequence and giant binaries and the 
integrated blue straggler sample. Here we include all binaries with 
P < 3000 days, and the last bin is normalized to account for the dif- 
ferent bin size. In the main-sequence plot, we also show the initial 
period distribution of the simulation within the same mass, period 
and spatial domain with the gray line. All simulated distributions 
are normalized by the respective NGC 188 sample sizes, and we 
show the true frequencies on the right y-axes. The respective NGC 
188 distributions are plotted in the thick-lined histograms. For 
comparison we also sh ow the Galactic field period distribution of 
Raghavan et al. (2010) with the dashed lines, normalized to the 
NGC 188 sample size using the field binary frequency within this 
period range. The simulation contains a large overabundance of 
short-period main-sequence binaries as compared to observations 
of both NGC 188 and the field. 

orbital parameters. Such changes in orbital parameters 
are accounted for, in the integrated sample, as separate 
systems with lifetimes that are equal to the duration at 
each orbital configuration. 

To close, Hurley et al. did not have access to RV data 
as comprehensive as ours for either the M67 or NGC 
188 binary population when constructing their simula- 
tion or analyzing their results. Therefore, we are able 
to investigate the accuracy of their model both in more 
detail and, perhaps more importantly, in observational 
dimensions not available to them. Many of the com- 
parisons discussed next truly test a priori predictions of 
their theory. 

9.2. Comparison of the Main-Sequence Hard-Binary 
Populations in NGC 188 and in the Simulation 

9.2.1. Frequency 

We find a MS solar-type hard-binary frequency at 4 
Gyr in the simulation of 64 ± 4 % (230/361) for P < 10 4 
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Fig. 13. — e — log(P) diagram for binaries in the IHurlev et all 
[|2005h simulation. We compare the main-sequence (top), giant 
(middle) and blue straggler (bottom) populations from the simula- 
tion. We show the 4 Gyr main-sequence and giant binaries. For the 
blue stragglers, we show the integrated (2-4 Gyr) sample, where 
each point represents a given blue straggler binary configuration 
during one time step. If a dynamical encounter alters the orbital 
period or eccentricity of a given blue straggler in a binary, the re- 
sult is shown by individual (often partially overlapping) points at 
each orbital configur ation, respectively. T he circularization period 
Pcirc = 15.0 days (Mathieu et al. 2004) is shown in the main- 
sequence e — log(P) plot by the dotted line. Note the unphysical 
"envelope" to the main-sequence plot that is not seen in the NGC 
188 ob servations (or those of oth er open clusters of a wide range 
in age. IMeibom fc M athieu 2005j). 

days. (Here the uncertainty is the expected counting er- 
ror given multiple realizations of the simulation.) We 
specify the binaries defining this frequency in the same 
manner as for the MS NGC 188 binaries in Section[3] The 
MS hard-binary frequency in the simulation has grown 
from its initial value of 47.9 ± 2.5 % (for the same mass, 
spatial and period domain). This increase in binary fre- 
quency is a result of the preferential evaporation of the 
solar-type single stars over their higher total-mass solar- 
type binary counterparts. 

The 4 Gyr MS solar-type binary frequency is twice that 
of the NGC 188 MS hard-binary frequency of 29 ± 3 %, 
and these two frequencies can be distinguished at the 
>99% confidence level. The simulated giant hard-binary 
frequency of 52 ± 16 % (11/21) is also higher than for 
the incompleteness-corrected NGC 188 giant hard-binary 
frequency of 31 ± 7 %, although given the large uncer- 
tainties these two values are not distinct at high confi- 
dence. These differences in binary frequency are directly 
linked to the high initial hard-binary frequency in the 
simulation, and the small subsequent change in frequency 
during the cluster lifetime. 



Fig. 14. — E ccentricity distribution for binaries in the 
IHurlev et al. (2005) simulation. We compare the main-sequence 
(top), giant (middle) and blue straggler (bottom) populations from 
the simulation, using the 4 Gyr main-sequence and giant binaries 
and the integrated blue straggler sample. Here we show only bina- 
ries with P c i rc < P < 3000 days, where t he circularization period 
Pc ir c = 15.0 days UMathieu et al.l I2004T) . In the main-sequence 
plot, we also show the initial eccentricity distribution of the simu- 
lation within the same mass, period and spatial domain with the 
gray line. All simulated distributions are normalized by the respec- 
tive NGC 188 sample sizes, and we show the true frequencies on the 
right y-axes. The respective NGC 188 distributio ns are plotted in 
the th ick-lined histograms. Finally, we show the Raghavan et al. 
(2010) distribution in the dashed lines, normalized to the NGC 
188 sample size using the field binary frequency within this period 
range. 

9.2.2. Orbital Period 

In Figure [T2"1 we plot the period distribution for the sim- 
ulation. We use the same gray-scale (color in the online 
version) coding for the simulated MS, giant and BS pop- 
ulations as we have done throughout the paper for NGC 
188. We normalize the simulated output by the respec- 
tive NGC 188 sample sizes and show the true frequencies 
on the right y-axes. For reference, we also include the 
NGC 188 distributions as the thick histograms and the 
normalized R10 field distributions as the dashed lines. 
As in Figure |6] we show binaries out to P — 3000 days in 
Figure [T2l and normalize the last bin to account for the 
different bin size. 

The MS period distribution of the simulation is signif- 
icantly different from that observed in NGC 188 (Fig- 
ure [6|). In particular, the simulation contains a large 
overabundance of short-period MS binaries as compared 
to both NGC 188 and the Galactic field. If we normalize 
the MS period distribution from the simulation to con- 
tain the same number of binaries as in the NGC 188 MS 
sample, a x 2 test shows that the two distributions are 
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Fig. 15. — Log period (top) and eccentricity (bottom) plotted against projected radius from the cluster center for the solar-type initial 
(left) and final (at 4 Gyr; right) binary populations. We show the main-sequence, giant and blue straggler binaries in light-gray circles, 
dark-gray diamonds and black triangles (green, red and blue in the online version), respectively. (I n the initial plots al l binaries are on the 
main-sequence.) On the top x-axes we convert parsecs to core radii using a core radius of 0.64 pc ijHurlev et al][2b05fl . Despite becoming 
more centrally concentrated, the overall form of the period distribution does not change over 4 Gyr of evolution. On the other hand, we 
notice a population of high-eccentricity binaries in the core of the cluster at 4 Gyr that were not present in the initial population. Their 
higher eccentricities are the result of dynamical encounters. 



distinct at the >99% confidence level. 

The distribution of binary periods has not changed sig- 
nificantly over the 4 Gyr of simulated evolution. A K-S 
test relating the initial and final hard-binary period dis- 
tributions (Figure [12]) shows no distinction. Thus this 
discrepancy between the simulated and observed popu- 
lations is a result of the flat initial distribution in log(P) 
used in the simulation. 

9.2.3. Orbital Eccentricity 

In Figure [T5| we show the e — log(P) diagram for the 
iHurlev eT al. (2005) simulation. The e — log(P) distribu- 
tion of the simulation MS binaries is strikingly different 
from the observed e — log(P) distribution of MS binaries 
in NGC 188 (top panel, Figure [3]), and indeed from ob- 
servations of many si milarly evolved open cluste rs and 
the Galactic field (see iMeibom fc Mathieul 120051 R10). 
The differences in the simulation include: 

1. The upper envelope on orbital eccentricity, rising 
with increasing period. This derives directly from 
the thermal initial conditions and the Kroupa tidal 
evolution prescription. The observed distribution 
of orbital eccentricities longer than the circulariza- 
tion period is better matched with a wide Gaussian 
in eccentricit y at all periods , as found in the young 
cluster M35 (Mathieij |2008) and in the field (R10), 
and as used in the modeling of detac hed binary 
evolution bv lMeibom fc Mathie"ul (|2005l ). 



2. The large number of circular binaries with periods 
greater than the observed circularization period of 
15.0 days. Many of these systems had circular or- 
bits initially due to the Kroupa pre-evolution pro- 
cess. An additional small population of long-period 
circular binaries was created through unstable mass 
transfer (common-envelope) events involving ini- 
tially non-circular MS binaries. 

3. The large number of eccentric binaries with periods 
shorter than the circularization period. This sug- 
gests that the tida l energy dissipation rate is un- 
dere stimated in thelHurlev et alj ( |2005l ) model. In- 
deed |Belcixnikre£2Dl200^~s^iggest that the con- 
vective tidal damping should be increased by a fac- 
tor of 50 to match observations of M67. 

Similarly, in Figure [14] we show the eccentricity dis- 
tributions from the simulation. Here we use the same 
period limits of P C irc < P < 3000 days as in Section [5] 
to define our sample, where P C irc = 15.0 days (for NGC 
188 found bv lMathieu et al.l 120041) . The marked excess 
in the number of circular orbits compared to NGC 188 
(Figured]) and the field (R10) is evident. This same ex- 
cess i s seen in comp arison with the young open cluster 
M35 (]Mathieull2008f ). Most of these circular orbits are in 
short-period binaries, and as such play a critical role in 
the formation of BSs in the simulation. 
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Fig. 16. — Secondary-mass (left) and mass-ratio (right) distributions for solar-type main-sequence binaries in the Hurley ct al. (2005) 
simulation (filled histograms) and in NGC 188 (thick-lined black histograms), both for binaries with P < 3000 days. We also show the 
initial populations from the model within the same mass and period range in the gray lines. We normalize the simulated distributions by 
the NGC 188 sample size, and show the true frequencies on the right y-axes. As the simulation only included stars with masses >0.1 Mq in 
the initial population, we show the secondary-mass distribution from 0.10 - 1.15 Mq and the mass-ratio distribution from 0.1 - 1. 



Finally, in Figure[T5]we show the initial and final distri- 
butions of eccentricity as functions of radius. Over the 
course of the cluster evolution, a population of higher- 
eccentricity binaries developed within ^5 core radii from 
the cluster center. (This population can also be seen 
in the e — log(P) diagram of Figure [131) No similarly 
high eccentricity binaries are found in the cluster halo. 
The increases in the orbital eccentricities of these bi- 
naries can be traced to dynamical encounters. These 
high-eccentricity binaries in the model are likely ana- 
logues of the high-eccentricity binaries observed in the 
core of NGC 188 (discussed in Section suggestive of 
a dynamical origin for some or all of the NGC 188 high- 
eccentricity binaries. 

9.2.4. Secondary-Mass and Mass-Ratio Distributions 

Unlike our observed NGC 188 sample where we re- 
quire a statistical method to derive the secondary-mass 
and mass-ratio distributions, within the simulation we 
know the secondary masses and mass ratios for every bi- 
nary. In Figure [TB] we show the 4 Gyr secondary-mass 
(left) and mass-ratio (right) distributions for the simu- 
lated MS binaries (filled histograms), and compare to 
the respective NGC 188 MS distributions (thick-lined 
black histograms) as well as the initial distributions in 
the simulation (thin- lined gray histograms). The excess 
in hard-binary frequency - both initial and at 4 Gyr - 
in the simulation compared to NGC 188 is again clearly 
evident. However here we focus on the shapes of the 
distributions. 

Both the secondary-mass and mass-ratio distributions 
of the simulation are essentially flat at 4 Gyr, except for 
the highest mass secondaries which are few in the sim- 
ulation. This flat shape is in marked contrast to that 
found in NGC 188, which rises to lower masses and mass 
ratios. The NGC 188 distribution also shows some ev- 
idence of a peak at comparable primary and secondary 
masses (" twins" ) , which is not seen in the simulation. 

In fact, the initial secondary- mass and mass-ratio dis- 
tributions of the simulation do rise to lower secondary 
mass and look much like the current NGC 188 distribu- 
tions. However, by 4 Gyr they flatten with subsequent 
dynamical evolution of the cluster. This would suggest 



that the NGC 188 initial distribution may itself have 
been even more peaked to lower masses than currently 
observed. However the simulation is silent with respect 
to the creation or evolution of a population of twins. 

The 4 Gyr mass-ratio distribution from the simulation 
is quite similar to that of the field binaries of R10, in 
that the distribution is consistent with being uniform 
over mass ratios of ^0.2 - 0.95. However the simulation 
does not reproduce the abundance of twins (q ~ 1), again 
suggesting that this feature seen in both the field and 
open clusters is likely linked to the processes of binary 
formation rather than dynamics. 

9.2.5. Summary 

The MS hard-binary population predicted by the sim- 
ulation is significantly different from the MS hard-binary 
population of NGC 188. Along some dimensions this re- 
flects on the adopted initial conditions of the simulation 
and along others it reflects on the physics of the simula- 
tion itself. 

The initial binary frequency and period distribution 
of the simulation are largely maintained throughout the 
cluster lifetime, and so the marked differences with the 
NGC 188 frequency and period distribution are the re- 
sult of inaccurate initial conditions. Specifically, the flat 
initial period distribution is distinct from the log normal 
period distribution of the observ ed young (~150 Myr) 
open cluster M35 (|Mathieull2008t ). and the initial hard- 
binary frequency integrated over period is a factor two 
higher than found in such clusters. Given that the initial 
binary period distribution imprints itself on the cluster 
throughout its evolution, it is imperative that the prop- 
erties of the initial binaries used in future TV-body simu- 
lations be guided by observations of such young clusters. 

The differences in the simulated and observed orbital 
eccentricity and e — log(P) distributions are also in large 
part the result of inaccurate initial conditions. Again, the 
eccen tricity distribut ions in o bserved open clusters (e.g. , 
M35. iMathieul 120081 : and see iMeibom fc Mathieul [2005ft 
argue that a Gaussian-like distribution in eccentricity, 
with a weak dependence on period above the initial cir- 
cularization cutoff, is a better choice than a thermal dis- 
tribution. 
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In contrast, the large frequency of eccentric binaries 
with periods shorter than the circularization period of 
M67 and NGC 188 indicates that the physics of tidal 
circularization incorporated in the simulation is not ac- 
curate. Arguably the tidal energy dissipation rate is too 
small, as suggested by lBelczvnski et aLl ((2008). This can 
be empirically tested with future simulations. On the 
other hand, the Kroupa (1995) tidal evolution prescrip- 
tion applied to the initial binary population evidently 
produced too many initial circular orbits at longer peri- 
ods. 

The presence in the simulation of a central concen- 
tration of binaries and especially binaries in the core 
with dynamically produced high orbital eccentricities is 
in agreement with the binary population of NGC 188 and 
is encouraging with respect to the dynamical encounter 
physics in the simulation. We have not attempted a spe- 
cific comparison of rates and frequencies, in part because 
iHurlev et al.l (2005) derive a core radius for the simula- 
tion of 0.64 pc at 4 Gyr, sm aller than observed in NGC 
188 ( 1.3 pc. lBonattoet all 120051 2.1 pc, IChumak et "all 
2010). This higher central concentration will change evo- 
lution timescales in ways that likely are not easily com- 
pensated for accurately. 

The correct initial secondary-mass and mass-ratio dis- 
tributions remain poorly known, and comparison with 
NGC 188 provides only a rather distant lever on the 
question. Based on the simulation and NGC 188 dis- 
tributions, it would appear that a secondary distribution 
rising to smaller masses is appropriate, and perhaps more 
so than used by Hurley et al. (2005). More intriguing is 
whether twins are present initially, given the hint of an 
excess of them existing both in NGC 188 and the field. 

The CMD was the primary interface with observation 
used by Hurley et al. (2005), and they were very success- 
ful in matching the observed CMD of M67. In particular 
they reproduced well the BS population of M67. Im- 
portantly, the BS population of the simulation derived 
largely from mergers of close binaries, which we find to 
be over-represented in the simulation due to the adopted 
initial binary population. A plausible hypothesis is that 
reducing the binary frequency by a factor of two and 
reducing the frequency of the closest binaries relative 
to those of longer period would substantially reduce the 
BS production rate and un der-predict the BS p opulation 
compared to M67. Indeed. IHurlev et~ al. (2005) note that 
their A-body models using initial perio d distributions 
such as lDuauennov fc Mayor] (|199lD and lKroupal (|1995f ) 
produce significantly fewer BSs than their M67 model 
(which used a flat initial period distribution). Thus such 
models may require an offsetting change in the physics of 
BS formation or the relative frequency of BS formation 
channels in order to recover agreement with M67. We 
begin discussion of this issue in the next section. 

9.3. Comparison of the Blue Straggler Populations in 
NGC 188 and in the Simulation 

9.3.1. Comparison of the Blue Stragglers in Binaries 

Using the integrated sample of BSs, we derive a BS 
hard-binary (P < 10 4 days) frequency of 29 ± 2 % 
(168/584). This hard-binary frequency for the BSs in 
the simulation is significantly lower than the NGC 188 
BS hard-binary frequency of 76 ± 19 %. In Figure[l7l we 
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Fig. 17. — Blue straggler hard-binary (P < 10 4 days) frequency 
as a function of time in the Hurley et al. (2005) Af-body simula- 
tion. To reduce the uncertainties and scatter we average every 
four adjacent time steps, so that each point represents ~240 Myr 
in time. We show the region in time covered by each point in the 
gray horizontal bars. For comparison, we also show the NGC 188 
BS binary frequency of 76 ± 19 % with the light gray band. 



plot the lHurlev et al.l (|2005[ ) BS hard-binary frequency as 
a function of cluster age (no longer limited to ages >2 
Gyr) . Over most of the lifetime of the simulated cluster 
the BS hard-binary frequency is less than half that of 
NGC 188 (and M67). Thus, although the simulation has 
an overabundance of hard MS binaries, they are decid- 
edly lacking amongst the simulated BS population. 

We estimate from the CMD location that the most 
massive NGC 188 BS in a binary with an orbital solu- 
tion has a mass of ~1.5 times the MS turnoff mass of 1.1 
M Q . The MS turnoff mass in the simulation is 1.7 M at 
2 Gyr and decreases to 1.3 M by 4 Gyr. Roughly 44% 
of the BSs in binaries (P < 3000 days) in the simulation 
(between 2-4 Gyr) have masses >1.5 times the respec- 
tive simulation turnoff mass. Thus, on average the BSs 
in the simulation are more massive than those in NGC 
188. 

In Figure fT8l we compare the NGC 188 BS secondary- 
mass and mass-ra tio distributions to those in the 
IHurlev et all (|2005h simulation. The BSs in the simu- 
lation do not show the strong peak at ~ 0.5 M Q that 
is observed in the NGC 188 BSs. Instead the distribu- 
tion is roughly uniform for secondary masses >0.2 M . 
Interestingly the observed and simulated mass-ratio dis- 
tributions for the BSs appear to be quite similar in form. 
This agreement between the mass-ratio distributions is 
the result of dividing high-mass companions by the high- 
mass BSs, both of which individually are inconsistent 
with the NGC 188 BSs. 

The binaries in the integrated BS sample of the sim- 
ulation also have significantly different periods and ec- 
centricities than the BS binaries of NGC 188. Specifi- 
cally, the simulation produces BS binaries with shorter 
periods and a higher frequency of large eccentricities. A 
2D K-S test comparing the observed and simulated BS 
e — log(P) distributions shows that the two distributions 
are drawn from different parent populations at the >99% 
confidence level. Likewise ID K-S tests comparing the 
observed and simulated eccentricity and period distribu- 
tions return distinctions at the 99% and >99% confidence 
levels, respectively. 

Thus the BS sample from the simulation has signifi- 
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Fig. 18. — Secondary-mass (left) and mass-ratio (right) distributions comparing the blue stra ggler binaries in the Hurley ct al. (2005) 
simulation to those in NGC 188. The thick-lined histograms reproduce the results shown in Figure [T0l for the NGC 188 blue stragglers. For 
the simulation, we use the integrated sample of all simulated BSs from 2-4 Gyr. The filled histograms show the simulated distributions 
normalize by the NGC 188 sample size. The true frequencies are shown on the right y-axes. For comparison of distribution shapes we also 
re-normalize the distributions from the simulation to contain the same number of blue straggler binaries as in NGC 188, shown with the 
dashed lines. 



cantly different binary properties from those of the NGC 
188 BS population. These discrepancies are almost cer- 
tainly linked, at least in part, to the inaccurate initial 
binary population used in the simulation. 

9.3.2. Comparison of the Blue Straggler Spatial 
Distributions 

The spatial distribution of the BSs in the simulation 
is also quite different from what we observe in NGC 188 
(Figure [TOl and see also Figure 9 from Paper 1). In Pa- 
per 1 we show that the NGC 188 BSs have a seemingly 
bimodal spatial distribution, with a population of 14 cen- 
trally concentrated BSs in the core and 7 BSs in the halo. 
The full NGC 188 BS population is not significantly more 
centrally concentrated than the solar-type single stars in 
the cluster. In the simulation, we do not observe this rel- 
atively large population of BSs towards the cluster halo. 
Instead the full BS population is centrally concentrated 
with respect to the single MS stars (at the 98% confi- 
dence level; also see Figure fl~5|) . Furthermore, if we only 
analyze the 4 Gyr BS population, all but one BS is within 
4.5 pc from the cluster center; the remaining BS is at a 
radius of 19.38 pc which is beyond the spatial extent of 
our NGC 188 survey. 

Other than the difference in number of BSs, the core 
and halo BS populations in NGC 188 are statistically 
indistinguishable. Both have consistent binary frequen- 
cies and indistinguishable distributions of orbital eccen- 
tricities, periods and mass functions. Furthermore both 
the core and halo populations contain one of the two 
short-period SB2 BSs, resp ectively. The distributio ns of 
vsini rotational velocities (|Mathieu fc G eller 2009|) are 
also consistent between the core and halo populations. 

Similar bimodal BS spatial distrib utions have also been 
observed in globular clusters (e.g.. iFerraro et al.l 1997 ; 
Zaggia et a l. 1997; Ferraro et alJ 2004: Sa bbi et al.ll2004 : 



r (core radii) 
9.4 12.5 15.6 



Lanzoni et a 



2007h . One interpretation for the origins of 



these peculiar spatial distributions has been the presence 
of two distinct formation mechanisms, with mass transfer 
domi nating in the halos a nd collisions dominating in the 
core (|Mapelli et al.ll2006r ). Alternatively the distribution 
may be a dynamical signature where encounters kick BSs 




Blue Straggler Stars 

Single Main Sequence Stars _ 

Single Giant Stars 



12 



14 



16 



r(pc) 



Fig. 19. — Cumu lative spatial distribution of stellar populations 
in the IHurlev et al.l (|2005D TV-body model. The blue stragglers 
from the integrated sample are shown in the solid line. This BS 
spatial distribution is compared to the single main-sequence and 
giant stars, shown in the dashed and dotted lines, respectively. 
The blue stragglers are more centrally concentrated than the single 
main-sequence stars (at the 98% confidence level). 



out of the core leaving them to either relax back towards 
the center or re main in the halo depen ding on the local 
relaxation time (jSigurdsson et al]|1994D . 

The similarities in the observed properties between the 
core and halo NGC 188 BS samples suggest that in NGC 
188 this spatial distribution is not the result of two dis- 
tinct formation mechan isms. Instead, as discussed in 
iGeller fc MaThleul ^OlH ). the majority of the NGC 188 
BSs are consistent with an origin in mass transfer. Dy- 
na mical encounters in the core similar to those suggested 
bv lSigurdsson et all (|1994l) (though not resulting in colli- 
sions) may be responsible for the BS spatial distribution. 
Alternatively the majority of NGC 188 BSs may have 
formed with minimal dynamical disturbance following a 
very similar spatial distribution to the binaries in the 
cluster. The "gap" now seen in the spatial distribution 
between radii of ~5 and ~15 arcmin may then be a result 
of mass segregation processes, and may mark the region 
where the relaxation time approaches the age of the BSs. 
(A similar m echanism has also be en s uggested by other 
authors, e.g., iMapelli et aT1l2004l and iLeigh eT~aUl2011l 
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Fig. 20. — Cumulative orbital period distributions comparing the NGC 188 blue stragglers to those in the Hurley et al. (2005) simulation. 
The NGC 188 blue straggler orbital periods are shown in the black points, with the SB2s circled. The period distributions from the simulation 
are divided by formation channel, with mass transfer in the solid line, mergers in the dotted line and collisions in the dashed line. In the 
left panel we show only blue straggler binaries in the integrated sample from the model with P < 3000 days, while on the right we no 
longer limit by period. In the right panel the gray-filled region marks periods beyond our observational completeness limit. The simulation 
predicts that merger products will either be found in contact binaries (in the process of merging) or in wide binaries that result from 
exchange encounters, both of which are inconsistent with the NGC 188 blue straggler periods. Collision products that remain in binaries 
are predicted to have a significantly lower frequency of binaries with P < 3000 days than observed for the NGC 188 blue stragglers. The 
mass-transfer products, on the other hand, all have periods within our completeness limit. 



to explain the bimodal BS spatial distributions observed 
in some globular clusters.) 

9.3.3. Testing Blue Straggler Formation Channels in the 
Cluster Simulation 

There are t hree m echanisms for forming BSs in the 
iHurlev et al.l (|2005T ) TV-body model: mergers, mass- 
transfer processes and stellar collisions. Here we inves- 
tigate the binary populations for BSs formed from each 
mechanism as they appear in the integrated sample, and 
compare their binary characteristics to those observed for 
the NGC 188 BSs. In addition to completed mergers, we 
include in the merger category any binaries that are cur- 
rently in contact or u ndergoing mass tra nsfer and even- 
tually merge. (In the lHurlev et al.ll2005l model, mergers 
result from Case A and/or B mass transfer. The model 
never contains a large enough population of triple stars 
to test BS formatio n through the recent K ozai-based 
merger mechanism of lPerets fc Fa brvckv 2009.) Between 
2 and 4 Gyr in the simulation, 75% (439/584) of the BSs 
in the integrated sample formed through mergers, 20% 
(117/584) of the BSs formed through collisions and 5% 
(28/584) of the BSs formed through mass transfer (that 
did not result in a merger). 

The relative rates of BSs created through each of these 
processes in the simulation are closely linked to the ini- 
tial binary population. The dominant formation channel 
for the BSs in the simulation is mergers, a direct con- 
sequence of the over-abundance of short-period binaries 
included in the primordial population. We've shown that 
the initial binary population used in the simulation is not 
appropriate for NGC 188; thus the absolute and relative 
BS formation rates cannot be expected to agree with the 
NGC 188 BSs. 

However the predictions for the BS binary properties 
of each formation channel (e.g. distributions of binary 
orbital parameters) are only minimally sensitive to the 
initial binary population, and therefore still provide a 
valuable means to investig ate the origins of t he NGC 
188 BSs. Here we use the IHurlev et ail (|2005t ) models 



to define theoretical expectations for properties of the 
BSs formed by each channel and compare to the ob- 
served properties of the NGC 188 BSs. Our goal is to 
use these properties to begin to determine which mech- 
anisms played significant roles in the formation of the 
NGC 188 BSs. 

We begin first with the BS binary frequency. As ex- 
pected, the merger mechanism predicts a low BS binary 
frequency (P < 10 4 days) of only 13 ± 2 % (61/439), 
much lower than the observed NGC 188 BS binary fre- 
quency. Furthermore, the majority of these binaries are 
contact systems, which we do not observe among the 
NGC 188 BSs. 

The collision mechanism predicts a total BS binary fre- 
quency of 74 ± 8 % (86/117). However, the vast ma- 
jority (74%; 64/86) of the BSs in binaries produced by 
collisions have periods longer than 3000 days, our com- 
pleteness limit for orbital solutions in NGC 188 (see Sec- 
tion HJ. 67% of the NGC 188 BSs have periods shorter 
than 3000 days, which contradicts the simulation result 
of only 19% (22/117) binaries with periods of less than 
3000 days from the collisional formation mechanism. 

The mass-transfer mechanism predicts a 100% BS bi- 
nary frequency, all of which are detached and have peri- 
ods of less than 3000 days. 

Next we consider the distributions of orbital parame- 
ters predicted by the simulation for binaries formed by 
each mechanism, beginning with the period distributions. 
The predicted cumulative period distributions for the 
three mechanism are shown in Figure 120} the left panel is 
limited to BS binaries with periods below our complete- 
ness limit while the right panel shows all BS binaries 
from the simulation. 

Two populations of merger products are evident in Fig- 
ure [201 a grouping of contact systems at short periods 
and a grouping at longer periods where the BSs were sub- 
sequently exchanged into binaries. The period distribu- 
tions of both of these merger populations are inconsistent 
with the period distribution of the NGC 188 BSs. 

The distribution of collisionally produced BS binaries 
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Fig. 21. — Cumulative orbital eccentricity d istributions compar - 
ing the NGC 188 blue stragglers to those in the lHurlev et al.l (120051 ) 
simulation, for binaries with P < 3000 days. The NGC 188 blue 
straggler orbital eccentricities are shown in the black points, with 
the SB2s circled. The eccentricity distributions from the simulation 
are divided by formation channel, with mass transfer in the solid 
line, mergers in the dotted line and collisions in the dashed line. 
The JV-body model predicts that collision products will have sig- 
nificantly higher eccentricities than observed in the NGC 188 blue 
straggler binaries. Both mass-transfer and merger products are 
predicted to have nearly all circular orbits. (The few non-circular 
orbits in these samples are the result of dynamical encounters.) 
The assumption of circular orbits for contact binaries is li kely valid. 
However the assumptions used in thejHurlcy ct al. (2005j) model of 
rapid circularization during mass transfer may no t be valid in all 
cases (Sokcr 2000; Bonacic Marinovic ct al. 2008; Scpinsky ct al. 
2009). Thus the accuracy of the predicted eccentricity distribution 
for mass-transfer products from the model is uncertain. 

with periods of less than 3000 days is consistent with 
the shape of the observed distribution of NGC 188 BS 
binaries. However, as already noted, the frequency of 
such binaries in the simulation is m uch too small. 

All of the mass-transfer BSs in the lHurlev et al.1 (=2005) 
model have periods within our detection limit, and are 
split between one population near 100 days (a result of 
Case B mass transfer from a red giant) and a second pop- 
ulation near 1000 days (a result of Case C mass transfer 
from an asymptotic giant). 7 Similarly all but one of the 
NGC 188 SB1 BSs are found with periods near 1000 days; 
the remaining SB1 BS has a period of ^120 days. The 
theoretical and observed distributions are not formally 
distinguishable. 

The NGC 188 SB2 BS binaries both have very short 
periods around 5 days. Interestingly, the only formation 
mechanism in the M67 model that produces BSs with 
such short periods is the merger mechanism. However 
both NGC 188 BSs are in detached binaries. Thus the 
M67 model struggles to produce BSs with similar char- 
acteristics to the two SB2 BS binaries in NGC 188. 

The observed distribution of orbital eccentricities is a 
challenge for all of the formation mechanisms. Collision 
products that retain a binary companion are predicted 
to have significantly higher eccentricities than observed 
for the NGC 188 BSs; the simulation result is rejected at 
the >99% confidence level (Figure [2Tj) . Conversely, the 
mass-transfer and merger mechanisms both predict that 
essentially all of the BSs in binaries should have circular 

7 T hough Figure 1201 appears to suggest that the IHurlev et al.l 
(2005) model favors Case B over Case C mass transfer, the differ- 
ence in frequency is instead due to two longer lived Case B products 
rather than a significantly larger number of BSs produced through 
Case B mass transfer. 



Fig. 22. — Cumulative blue straggler mass distr ibution com- 
paring the NGC 188 blue stragglers to those in the IHurlev et al.l 
(2005) simulation. We normalize all masses by the respective main- 
sequence turnoff mass to facilitate the comparison. The integrated 
sample of blue stragglers from the simulation is limited here to 
only include binaries with P < 3000 days to reflect our observa- 
tional completeness in orbital solutions. Our NGC 188 blue strag- 
gler mass estimates are shown in the black points, with the SB2s 
circled. The primary mass distributions from the simulation are 
divided by formation channel, with mass transfer in the solid line, 
mergers in the dotted line and collisions in the dashed line. The 
simulation predicts that collision products are the most massive, 
followed by mergers and then mass-transfer products. Both the 
collision and merger products are significantly more massive than 
the NGC 188 blue stragglers, while products of mass-transfer have 
a mass distribution that is statistically indistinguishable from that 
of the NGC 188 blue stragglers. 

orbits. In the model, any non-zero orbital eccentricity 
prior to the BS formation is quickly circularized by tides 
during the mass-transfer or contact binary stages. The 
few non-circular orbits in both of these samples are the 
results of subsequent dynamical encounters. 

Three of the SB1 BSs in NGC 188 do have orbits con- 
sistent with being circular. Theoretical main-sequence 
tidal circularization timescales for binaries at ~1000-day 
periods are pr ohibitively long and canno t explain these 
circular orbits ([Meibom fc Mathieul l2~005') . Indeed, there 
are no circular MS binaries with periods near 1000 days 
in NGC 188. Thus these BS circular orbits may be evi- 
dence of a mass transfer origin. 

However, the other 10 long-period BS binaries in 
NGC 188 have measurable eccentricities, one as high 
as 0.8. Recently, observational and theoretical evi- 
dence has begun to suggest that mass transfer, partic- 
ularly involving an asymptotic-giant star, will not al- 
ways lead to rapid circularization. Examples of ec- 
centric post-m ass-transfer systems may include barium- 
star systems (|Jorissen et al.l 119981) , post- AGB binaries 
(jvan Belle et al.Hl999L and references therein), and pos- 
sibly also the long-per i od ec centric BSs in the galac- 
tic halo ([Carney et al.l 120051 ). Theoretical work has 
suggested mech anisms that can yield "eccentricity- 
pump ing (e.g. ISokerl l2000l : 1 onacic Marinovic et al.l 
|2008^ J Scpinsk y _et al.ll2009f) . which were not included in 
the IHurlev et al.l (|2005f) model. Thus at this time the 
implications of the observed eccentricity distribution for 
the mass-transfer channel remains uncertain. 

The eccentricities of the two short-period SB2's are 
interesting. The double BS binary 7782 found in the 
cluster halo has a circular orbit, while the BS-MS pair 
5078 has a small eccentricity of 0-1 7 . Bot h are discussed 
in more detail in lMathieu fc Gelled (j2009ft . 7782 is likely 
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the result of a dynamical encounter and exchange that 
may have ejected it to the halo; presumably after such 
an encounter it would have at least initially had a highly 
eccentric orbit. On the other hand, 5078 has managed to 
sustain a small eccentricity against a theoretically short 
tidal circularization timescale. 

In the previous section we found that the BSs in bi- 
naries within the simulation are generally more massive 
than those in NGC 188, where all are <1.5 times the 
turnoff mass. If we divide the integrated BS binary sam- 
ple from the simulation by formation mechanism, we find 
that the mean masses from mergers, collisions and mass 
transfer are 1.43 ± 0.05, 1.72 ± 0.10 and 1.3 ± 0.03 times 
the turnoff mass, respectively. Collision products in bi- 
naries are predicted to be the most massive BSs, followed 
by merger products and finally mass-transfer products. 

In Figure [22] we compare the mass distributions from 
the simulated BSs to that of the NGC 188 BS binaries 
(with orbital solutions). K-S tests show that the NGC 
188 BSs are significantly less massive than predicted by 
the collisional and merger hypotheses, at the >99% and 
95% confidence level s, respectively. T h is find ing is con- 
sistent with that of IGlebbeek fc Pols! (|2008h . who rule 
out recent formation of massive blue stragglers through 
collisions in NGC 188. Thus both collisions and mergers 
generally produce more massive BSs than those observed 
in NGC 188. On the other hand, a similar comparison 
between the predicted BS masses for mass-transfer prod- 
ucts to the NGC 188 BS masses shows no significant dis- 
tinction. 

Finally, and f or co mpleteness, we note that 
iGeller fc Mathieul (|2011l) present an analysis of the 
secondary-mass distributions of the SB1 BS binaries of 
NGC 188 (e.g., Figure fT0|) and those formed in a new 
TV-body simulation of NGC 188, created using a very 
similar code to that used by iHurlev et al.l (|2005jl , but 
with a more realistic initial binary population. They find 
that the observed secondary-mass distribution, narrowly 
peaked around a mean mass of 0.53 Mq, is not consistent 
at >99% confidence with the theoretical secondary-mass 
distribution of BS binaries formed by collisions, which 
have a mean mass nearly twice that observed. On the 
other hand, the observed secondary-mass distribution 
is fully consistent with a population of carbon-oxygen 
WDs with masses between 0.5 M Q and 0.6 M©. They 
also find that formation of BSs in hierarchical triples 
via the Kozai-induced merger mechanism is not favored 
by the observed distribution, but is only rejected at the 
98% confidence level. 

9.3.4. Summary 

To summarize, our comparisons show that BSs formed 
by collisions that retain a companion will have a signif- 
icantly lower binary frequency (P < 3000 days), longer 
orbital periods, higher orbital eccentricities, higher 
masses and higher-mass companions than are observed 
for the NGC 188 S B1 BS binarie s. BSs formed by merg- 
ers as modeled by IHurlev et al.1 (|2005[ ) are also highly 
inconsistent with the observations, most notably having 
a significantly lower binary frequency, shorter periods, 
and higher masses than are observed. Thus we can con- 
fidently rule out both mergers within isolated close bina- 
ries and collisions as the dominant formation channel for 
the NGC 188 BS binaries. Predictions for mass-transfer 



products are consistent with most of the binary proper- 
ties of the NGC 188 SB1 BSs, which comprise two-thirds 
of the NGC 188 BSs. 

We are unable to inve stigate the Ko z ai-ind uced merger 
mechanism within the IHurlev et al.l (|2005t l model, as 
this model never contains sufficient triples. Importantly 
the Kozai-induced merger and mass transfer hypothe- 
ses make additional observable predictions that can dis- 
tinguish between the two mechanisms. Mass trans- 
fer predicts WD companions while the Kozai-induced 
merger mechanism predicts predominantly MS compan- 
ions. Relatively young (< 0.4 Gyr) WD companions 
are hot enough to contribute a significant excess in UV 
emission to the combined spectrum of the BS binary, 
and can be detected through HST observations. Ad- 
ditionally mass-transfer products are expected to show 
a higher d epletion of lithium th a n merger and collision 
products (iLombardi et all 119951: ISandauist et all 119971; 
iNorris et al.l Il997t IGlebbeek et al.l 12010^ (The field 
BS binaries for which lithi um abundances are available 
all show lithium depletion; Ca rney et al.ll200~il ) Mass- 
transfer processes^ia^aJso_r^ducemodified CNO abun- 
dances (jShetrone fc Sandauistl 120001) and enhancements 
in s-process elements ([Carney et al.ll2001f) from the pro- 
cessed material of the giant star donor. Such observa- 
tions of the NGC 188 BS population are currently under- 
way and will soon help to further constrain the origins 
of the NGC 188 BSs. 

However mass transfer cannot explain the compan- 
ion masses or orbital periods of the NGC 188 SB2 
BSs. These two BS b inaries along with S0182 of M67 
(jSandauist et al.ir2003[ ) represent intriguing case-studies 
for future A-body simulations. We suspect that they 
are among the most direct observational evidence for the 
close dynamical encounters of binaries that TV-body stud- 
ies have predicted for decades. 

Additionally there are five BSs in NGC 188 that aren't 
detected as binaries. The M67 model suggests that these 
BSs likely did not all form through mass transfer, while 
both collisions yielding very wide companions and merg- 
ers predict such BSs. 

Finally, we anticipate that shifting the MS binary pe- 
riod distribution away from short-period binaries (to bet- 
ter match the observations of MS binaries) will shift the 
dominant formation mechanism away from mergers and 
thereby change the properties of the resultant BS pop- 
ulation, perhaps in better agreem ent with the binary 
observations. But recall that the IHurlev et al.l (|2005l ) 
simulation successfully matched the number and CMD 
distribution of the M67 BSs. To retain this agreement 
in BS frequency with a reduced number of initial short- 
period binaries may require modification of the physics 
and efficacy of the formation channels. We will explore 
this in a future paper discussing a more accurate A-body 
simulation of NGC 188. 

10. SUMMARY AND CONCLUSIONS 

Through a series o f papers (Papers 1 , 2 , 
iMathieu fc G^ne^l2009l and IGeller fc Mathieul I20H of 
which this is the most recent contribution, we study the 
dynamical state of the old (7 Gyr) open cluster NGC 
188. In previous papers we identify a complete sample of 
single and binary cluster members among the solar-type 
stars in NGC 188, present orbital solutions for the 
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majority of our detected binaries, and study the BSs in 
detail. Here we analyze the hard-binary population of 
the cluster, including the hard-binary frequencies and 
distributions of binary orbital parameters of the MS, 
giant and BS populations. We compare these results in 
detail to the observed Galactic field binary population 
of RIO as well as the simulated binary population in 
the lHurlev etail (|2005l ) iV-body open cluster simulation 
to study the impact that relaxation processes and close 
dynamical encounters may have on the cluster's binary 
population. 

We find the properties of the solar-type MS binaries 
in NGC 188 to be consistent with similar binaries in 
the Galactic field. The NGC 188 MS binaries have a 
log-period distribution that increases towards our detec- 
tion limit and is consistent with the log-normal distribu- 
tion found by RIO for solar-type binaries in the Galactic 
field. The MS eccentricity distribution has a roughly 
Gaussian form (for P C i rc < P < 3000 days), also con- 
sistent with RIO. Likewise the MS binary secondary- 
mass and mass-ratio distributions are indistinguishable 
from those of similar binaries in the Galactic field. 
In NGC 188 both the secondary-mass and mass-ratio 
distributions appear to rise to wards lower masses and 
mass ratios, as also observed bv lMazeh et al.l (|2003[ ) and 
iDuauennov fc Mayor! (j 19911 ) in the field. Both NGC 188 
distributions also show evidence for a secon dary peak 
at ma ss ra tios near unity, as o bserved by RIO. iTokovininl 
(|2000f ) and lFisher et~all (|2005l ) in the field. The observed 
mass-function distribution for the NGC 188 MS binaries 
is formally consistent with a secondary-mass distribu- 
tion where sec ondaries are drawn from either a standard 
IKroupal (l200l IMF or a uniform mass-ratio distribution. 

The observed global MS hard-binary frequency in NGC 
188 is 23 ± 2 %, which when corrected for our incom- 
pleteness results in a true MS hard-binary frequency of 
29 ± 3 % for P < 10 4 days. This hard-binary frequency 
is marginally higher than that of the RIO field within the 
same period range, of 19 ± 2 %. This somewhat higher 
hard-binary frequency in NGC 188 may be the result of 
dynamical relaxation processes where lower total mass 
single stars are preferentially ejecte d from the cluster, a s 
is seen in TV-body models such as iHurlev et al.1 ((2005). 
Two-body relaxation is evident in the central concentra- 
tion of binaries in the cluster. 

Additional evidence for dynamical encounters may be 
seen in the population of high-eccentricity binaries in the 
core of the cluster but not in the halo (Figure [5]) . High- 
eccentricity binarie s such as these are created within the 
IHurlev et al.l (|2005t ) model in the cluster core via dynam- 
ical interactions throughout the cluster lifetime. 

Importantly, the NGC 188 BS binaries are significantly 
different from the NGC 188 MS solar-type binaries, pro- 
viding some of the most distinguishing characteristics of 
the BS population found to date. The observed BS hard- 
binary frequency in NGC 188 is 76 ± 19 %, three times 
the observed solar-type MS hard-binary frequency. The 
BS and MS hard binary frequencies are distinguishable 
at the >99% confidence level. 

The BS binary e — log(P) distribution is distinct from 
that of the MS at the 99% confidence level. The BS SB1 
binaries in NGC 188 have particularly remarkable orbital 
characteristics. All but one have periods near 1000 days, 
most have modest eccentricities (with three in circular 



orbits), and their secondary- mass distribution is narrow 
and peaked, with a mean mass of ~0.5 Mq (Figure IT0|) . 

The two short-period (P < 10 days) BS binaries 
are both observed as SB2s. These SB2 BSs have the 
most massive companions to any of the BSs in NGC 
188. Indeed BS ID 7782 has a mass ratio of unity 
and is composed of two BSs. We argue that these 
short-period BS binaries likely have a dynamical origin 
(|Mathieu fc Gelled I20091) . 

SB2 BS, ID 5078, has a companio n star on the upper 
MS. Normal stellar evolution tracks (|Marigo et al.ll2008f) 
predict a mass for the BS primary star that is 15% to 
30% lower than the mass we derive using the kinematic 
mass ratio. Thus normal stellar evolution tracks may 
underestimate the masses of BSs. 

Com parisons to the BS population in the lHurlev et all 
(f2005h TV-body open cluster model show that we can rule 
out the hypothesis that the NGC 188 SB1 BS binaries 
have a collisional origin. Collision products that retain 
binary companions are predicted to have significantly 
higher eccentricities, higher masses, higher-mass com- 
panions, and a lower frequency of binaries with P < 3000 
days than observed for the NGC 188 SB1 BS binaries. 

Likewise mergers within isolated binaries will not re- 
produce the binary properties of the NGC 188 BSs. BSs 
produced by isolated mergers are too massive, have a bi- 
nary frequency that is much too low, and the few binaries 
have too short of periods. 

Predictions from the mass transfer mechanism are con- 
sistent with most of the binary properties of the NGC 188 
SB1 BSs, which comprise two-thirds of the NGC 188 BSs. 
Mass transfer products are predicted to have a very high 
binary frequency, with long orbital periods, and compan- 
ions of about 0.5 Mq, all close ly consistent with the NG C 
188 B Ss (as also discussed in lGeller fc Mathieu|[201ll ). 

The IHurlev et ail (|2005[ ) model never contains a suffi- 
cient frequency of triples to investigat e the Kozai-induced 
merg er mechanism for BS f ormat ion ([Perets fc Fabrvckvl 
120091) . iGeller fc Mathieu! (|2011[) find that the observa- 
tions do not favor this formation mechanism, but the 
data only reject the hypothesis at the 98% confidence 
level. Spectroscopic observations to determine the abun- 
dances of the SB1 BSs as well as HST observations aimed 
at detecting the FUV flux from the WDs predicted by the 
mass transfer mechanism are currently underway, and 
will further constrain the origins of the NGC 188 BSs. 

Comparing detailed observations of open cluster bi- 
naries with predictions from sophisticated TV-body open 
cluster models is a powerful method for investigating the 
dynamical evolution of a binary population and the for- 
mation of the BSs. Comparisons with observations, like 
those presented in this paper, are also essential to de- 
termine the validity of the models. In this paper we 
show that accurately defining the initial binary popula- 
tion is critical. We suggest that this is best accomplished 
through comparisons with detailed observations of binary 
populations in young open clusters, although every indi- 
cation to date is that the field binary population is a 
reasonable proxy. iV-body simulations modeling NGC 
188 and using an observationally defined initial binary 
population are underway and will soon allow for a de- 
tailed study of the dynamical evolution of the old open 
cluster NGC 188. 
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